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ABSTRACT 

-In  general,  the  investigation  of  practical  dynamical  systems  is  more 
heuristic  than  rigorous;  consequently,  most  literature  attempts  to 
provide  a conceptual  key  to  the  general  treatment  of  the  response  by 
correlation  with  linear  theory.  Practically  none  of  this  literature 
sets  forth  any  new  ideas  of  basic  methods  of  attack.  The  usual  attack 
is  by  extending  present  classical  methods  rather  than  by  inventiyig  new 
basic  approaches.  The  major  limitation  found  in  the  existing  work  is 
the  lack  of  comprehensive  understanding  of  the  basic  parametei's  of 
simple  nonlinear  oscillators. 

This  thesis  presents  an  accurate  solution  of  several  types  of  damping 
on  the  dynamical  behavior  of  harmonically  forced  single-degree-of- 
freedom  systems.  The  study  is  based  on  the  fundamental  parameters  of 
the  system.  The  fundamental  response  is  described  with  reference  to 
the  equivalent  damping  energy,  the  Ritz  method,  and  dimensional  analysis. 
Dimensional  analysis  is  used  to  develop  a method  for  predicting  the  gen- 
eral response  diagram  characteristics.  The  high  accuracy  of  the  solutions 
permitted  the  collection  of  some  very  important  design  data.  The  results 
are  presented  in  both  graphical  and  tabular  form  and  may  be  useful  to  those 
engaged  in  calibration,  design,  and  data  analysis  of  work  requiring  an 
accurate  solution.  Also,  the  linearized  methods  give  results  sufficiently 
accurate  for  many  engineering  applications.  The  prediction  of  response 
characteristics  by  dimensional  analysis  should  be  of  interest. 

The  time-histories  of  displacement,  velocity,  and  acceleration  are  presented 
along  with  the  response  diagrams  of  displacement,  velocity,  acceleration, 
phase  angle,  and  energy  curves  by  both  the  linear  and  accurate  solutions. 
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The  following  nomenclature  is  used  in  this  thesis: 

A Frontal  drag  area,  ft“^ 

C General  damping  coefl’icient  3 , units  vary 
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This  thesis  is  an  analytical  and  n'jmericai  investigation  of  the 
effects  of  several  different  types  of  damping  on  a harmonically-forced  single- 
degree-of-freedora  system.  At  present  there  exists  no  basic  comparison  tech- 
nique or  analytic  proced'ure  for  the  accurate  plotting  of  the  response  of  a 
damped  physical  system.  The  main  objective  is  the  establishment  of  the 
physically  correct  expression  for  the  damping  ratio  so  that  the  effects  of 
the  damping  function  can  be  accurately  described  and  compared  with  other 
functions.  The  customary  critical  damping  ratio  does  not  describe  cases 
other  than  the  viscous  case.  This  thesis  presents  the  correct  ratios  for 
four  basic  damping  functions  and  an  approach  for  determining  the  e.xpresGb  n 
for  a general  damping  function. 

The  effects  of  damping  on  the  general  dynamical  response  of  a physi- 
cal system  is  inportant  for  both  qualitative  and  quantitative  investigations. 

Qualitative  studies  are  important  because  of  their  wide  application 
to  distributed  systems.  In  general,  they  are  treated  by  applying  assumptions 
that  will  simplify  the  problem  analytically.  The  conceptual  key  to  the  gen- 
eral idealization  of  distributed  systems  is  the  understanding  and  intuitive 
knowledge  of  the  problem.  Once  the  problem  is  well  defint-d  physically,  avail- 
able mathematical  tools  may  be  applied.  Almost  without  exception  all  vibi-aticn 
problems  are  treated  analytically  as  lumped-parameter  systems. 

The  analytic  procedures  describe  the  system  in  terms  of  linear  theory. 
This  description  generally  takes  the  form  of  successive  linear  approximations. 
The  correct  assumption  for  the  solution  almost  never  is  obvious.  For  example, 
in  the  dynamic  calibration  of  shock  and  vibration  pickups,  measured  outputs 
are  sometimes  differentiated  to  determine  velocily,  acceleration,  and, nr  jerk 
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of  a sinusoidal  motion  applied  to  a system.  The  assump’-ion  is  that  there 
exist  no  higher  harmonics  in  the  motion.  The  basic  difficulty  with  this 
differentiating  process  is  that  a high-frequency  component  can  have  an 
arbitrarily  small  amplitude  in  the  original  wave  form  and  an  arbitrarily 
large  amplitude  in  the  differentiated  wave  form,  if  the  frequency  of  the 
component  is  high.  The  component  frequencies  may  be  either  amplitude, 
phase,  and/or  frequency  distorted.  The  effect  of  differentiating  an 
assumed  linear  motion  to  obtain  the  solution  would  be  questionable.  AiS(  , 
these  steady-state  solutions  are  obtained  only  if  the  physical  paranieterc 
change  in  a specified  way. 

Almost  all  of  the  analytical  and  experimental  work  has  very  little 
basis  for  comparing  the  effects  of  damping  with  other  systems.  This  work 
usually  is  limited  to  a study  of  eunplitude,  velocity,  acceleration,  or 
phase  angle  instead  of  a simultaneous  evaluation  of  all  the  parameters. 
This  is  particularly  true  for  nonlinear  response.  Indeed,  some  authors 
neglect  the  effects  of  damping  or  give  only  a cursory  description  of  its 
effect  on  the  system. 

In  order  to  make  a complete  qualitative  and  quantit.al ive  analysis 
of  the  problem,  every  pertinent  physical  propei’ty  should  be  described  in 
detail.  The  general  and  specific  valuations  of  t.hese  prcpertles  must  be 
accur-ately  described. 

The  small  amount  of  experimental  data  available  is  It  r specific  ei, 
vironmentai  conditions  and  is  taken  willi  met. hods  and  equipment  that  an.  tu' 
well  defined.  Also,  it  is  difficult  to  determine  tne  effect  . f intierent 
equipment  characteristics.  For  example,  an  acre leromet  ei--cabl e comt  inati>. 
can  completely  filter  harmonic  componetils.  Tiit‘  usua.  methods  lead  to  an 
analysis  of  the  maxim\im  values  of  the  resp.’tise;  Uiei-.  l'  re,  full  i!..:ig!it  in 
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the  problem  is  not  obtained.  A detailed  analysis  of  the  motion  should  in- 
clude a study  of  the  time-histories  of  the  response.  In  order  to  accomplish 
this  it  is  necessary  to  have  an  extremely  accurate  solution  to  the  problem. 

Since  most  of  the  basic  theory  of  operation  of  physical  systems  rests 
on  the  assun^ition  that  the  linear  theory  with  constant  coefficient  describes 
them,  it  is  logical  to  evaluate  the  problem  in  this  manner. 

The  most  logical  attack  to  the  problem  would  be  to  obtain  the  exact 
steady- state  and  transient  solution.  The  most  practical  approach  would  be 
to  do  this  for  the  most  common  types  of  damping,  i.e..  Coulomb,  viscous, 
velocity- squared,  and  displacement- squared,  that  are  encountered  in  practical 
applications.  Since  a multiple  degree  of  freedom  system  would  itjtroduce 
other  unknowns,  a sinusoidally  forced  single-degree  of  freedom  syste.m  will  be 
considered.  This  will  provide  the  information  necessary  for  an  accui-ate 
evaluation  of  the  effects  of  damping  on  system  response. 


The  object  of  this  thesis,  therefore,  is  to  (l)  obtain  an  ai-.-iiratt- 
solution  to  the  four  basic  damping  cases,  (2)  compare  the  results  with  tlie 
basic  theory,  (3)  investigate  the  general  and  specific  effects  of  dainpinr. 
(U)  reach  conclusions  concerning  these  effects.  (5)  obtaiii  accurate  desigi. 
data,  and  (6)  provide  an  adequate  base  for  future  research  of  combitied  iami 
ing  cases. 


CliAPTER  I 


STATEMENT  AND  FORM’JLATION  OF  THE  THEORETICAL  PROBLEM 

The  dissipation  of  energy  of  some  type  and  degree  occurs  in  ail 
physical  systems.  The  effect  on  mechanical  systems  is  regarded  generally 
as  being  undesirable  except  for  vibration  control  or  shock  isolation  prob- 
lems. The  desirable  effects  of  damping  have  been  used  to  good  advantage 
in  some  limited  instrument  design  areas.  For  example,  the  mechanical  move- 
ment of  graphic  recorders  may  be  damped  in  such  a manner  as  to  preserve  the 
phase  shift  characteristics  over  a wide  frequency  range.  It  seems  possible 
that  other  desirable  characteristics  of  the  dissipation  forces  coaid  be 
utilized  if  their  characteristics  were  knowii. 

In  order  to  fully  utilize  the  op>timum  response  characteristics  of 
a dynamical  system  it  is  necessary  to  know  such  steady-state  properties  as 
(a)  resonant  bandwidth,  (b)  damping  ratio,  (c)  damped  natural  frequency, 

(d)  phase  angle,  (e)  phase  angle  shift,  (f)  harmonic  content  of  the  .motion. 

(g)  the  maximum  displacement,  velocity,  acceleration,  and  (h)  tlie  wavefoi-m 
of  the  response. 

Pi-actically  all  mechanical  apparatus  are  subjected  to  acceleral  1 - !i 
loading.  The  evaluation  of  the  transient  response  to  this  ioadi!ig  is  alm.a't 
entirely  experimental.  The  effects  of  shock  are  determined  by  subjo'-'t  ing 
the  equipment  to  a known  impact  condition  and  observing  the  failure  Ly  vi.sual 
inspection  or  photographic  techniques.  This  pi'ocediu’e  idiviously  is  is!desl2-- 
able  since  it  would  be  belter  to  bc>  able  to  predict  the  (,'ffocts  in  the  initial 
design  stage.  The  basic  difficulty  is  evaluating  the  effects  of  rate  of 
loading  on  the  test  piece  and  the  t.est  ifist  rsimentatlen.  Th<'  < va;  u.ai  i.  f 
the  rate  of  loading  for  various  t,\qu-s  of  drnnpitig  fuiictions  should  in-1.  !■ 

If 


W 


analysis  of  the  third  derivative  with  respect  to  time  (jerk).  A study  of  these 


data  could  provide  sufficient  information  to  design  suitable  test  instriumentation. 


In  shock  isolation  problems  it  is  desirable  to  coraj;are  the  response 


characteristics  of  a damped  system  to  determine  the  peak  a.mpditude,  velocity, 


and  acceleration  versus  time  of  occurrence.  The  simultaneous  evaluation  of 


transient  energy  would  be  required  to  choose  the  optimum  system.  Since  the 


wave  forms  of  displacement  and  velocity  vary  with  the  type  and  degree  of 


damping,  the  energy  will  vary  considerably. 


The  establishment  of  a suitable  system  of  equations  that  will  allow 


a qualitative  and  a quantitative  study  of  these  effects  requires  careful  con- 


sideration. 


The  qualitative  aspect  of  the  proble.m  requires  the  establishment  of 


parameters  that  can  be  applied  to  a physical  problem.  Sound  judgment  as  to 


how  far  an  idealization  may  be  carried  without  affecting  the  response  must 


be  established.  However,  actual  design  and  practical  applications  require  a 


complete  and  accurate  knowledge  of  the  general  response  of  the  system,  which 


is  obtainable  only  through  quantitative  analysis.  The  above  statements  indi- 


cate a rather  broad  attack  on  the  problem;  however,  practical  limitations  re- 


strict the  number  of  equations  considered.  Therefore,  these  investigatioru 


are  limited  to  those  damping  forces  which  are  basic  to  design  practice. 


These  damping  forces  are  CouJomb,  velocity,  velocity- squared,  and  displace- 


ment-squared damp)ing. 


A wide  range  of  problems  can  be  rep>resented  analytically  with  con- 


stant mass  and  constant  spring  piropjerties . This  is  p'articul  ally  lime  in  inst.r  .- 


merit  design  since  spu’ings  with  hysteresis  errors  of  no  more  t.han  .0.'  p r -en 


accui’acy  are  reaiiily  available.  Therefore,  constant,  mass  atm  sp  ritig  chat - 


acteristlcs  can  be  considi'red  without  loss  of  gr-neral  it.v 


oupierS'cripl.  nimibcrs  refer'  to  the  Ref.'t'encer . p.  plii. 


\ 
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Almost  none  of  the  forcing  functions  are  sinusoidal;  however,  many 
are  periodic.  The  effect  of  the  forcing  function  on  a general  study  of 
dissipation  effects  is  not  significant,  since  any  function  can  be  completely 
described  by  amplitude,  period,  and  wave  form.  Therefore,  a sinusoidal 
function  will  be  considered. 

The  system  to  be  studied  is  represented  schematically  in  Fig.  1. 

It  consists  of  a mass  M,  a linear  spring  of  rate  k,  a harmonic  forcing  futiction 
F sin(um),  and  a generalized  damper  represented  by  Cf(x,x).  Equating  the  su-m 
of  the  non- inertial  forces  in  the  direction  of  moticn  on  the  mass  to  mass 
times  acceleration  and  rearranging  yields 

Mx  + Cf(x,x)  + kx  = F sin(u)t)  . (l) 

In  Eq.  (l)  kx  is  the  spring  force  on  the  mass,  cf(x,x)  is  the 
damping  force,  Mx  is  the  inertia  force,  and  F sin(cut)  is  the  forcing 
function.  The  instantaneous  displacement  and  time  are  x and  t,  respec- 
tively. The  forcing  function  eimplitude  is  F and  uu  is  the  circular  fre- 
quency of  the  forcing  function. 

Excluding  the  cases  of  no  damping  and  viscous  damping,  there  exists 
no  exact  analytic  method  of  solving  the  above  equation  that  will  give  both 
steady-state  and  transient  solutions.  The  solutions  that  are  available 
apply  only  to  special  cases  and  then  the  linear  apprtixi.mations  are  to  "first 
order"  accuracy. 

Equation  (l)  approximates  many  real  systems;  therefore,  the  evalua- 
tion of  the  effects  of  basic  danqlng  functions  should  r*'veal  useful  design 
information.  Also,  a wide  variety  of  problems  are  predominaritiy  affected 
by  one  of  the  basic  forms  of  damjlng.  For  exarajle.  the  fluid-dyn.'imic  drag 
on  a fully- submerged  body  is  predominantly  velocity-squared. 


System  Diagram 


t X 

+ X 
-t  X 


Mx  + 9^  + kx  = F sin  u)t 

Mi?  + B-  X + kx  = F sin  cat 
Mx  + 0.^x|x|+kx-F  sin  cat 

J 

• * 2 

Mx  + 8j^  X -t  kx  - F sin  cat 


Equations  of  Motion 
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Veiocity- squared  damping  is  most  predominant  ir.  practical  applications 

.2 

of  fluid  dynamics.  The  magnitude  of  the  damping  force  is  C^pA  . 

For  a particular  problem  all  of  the  variables  are  well  defined 
except  C^.  The  value  of  can  be  obtained  by  testing  in  a wind  or  water 
tunnel  or  by  reference  to  previous  work.  Hoerner^  has  published  a book 
that  gives  the  values  of  for  many  configuration,  ncountered  in  prac- 
tice. The  combination  of  these  coefficients  with  the  solution  to  Eq.  (l) 
will  permit  an  accurate  description  of  the  motion  of  systems  affected  by 

7 

this  damping  force.  Also,  fluid  flow  in  pipes  has  been  represented  by  a 

.1.8  II 

constant  times  x . The  fundamental  study  of  a Helmholtz  resonator  is 

closely  represented  by  Eq.  (l)  except  for  the  radiation  losses  at  the  end 

of  the  neck. 

Coulomb  friction  is  a constant  force  and  independent  of  velocity 
or  displacement.  The  direction  of  the  force  opposes  the  relative  motion 
across  the  damper.  For  periodic  motion  the  relative  velocity  is  zero  twice 
during  each  cycle;  in  practice  the  coefficient  of  friction  is  a function  of 
the  materials  at  the  interface.  The  static  and  dynamic  coefficierits  differ 
by  a factor  of  2 to  5,  and  an  effective  value  of  the  coefficient  is  custom- 
arily used.  The  coefficient  of  rolling  friction  for  preloaded  tall  and 
roller  bearings  is  fairly  constant.  In  general,  shafts  sliding  Ihrougli 
sleeve  bearings  and  structural  join ts  imay  be  appro.ximatod  with  an  effc.-- 
tive  damping  coefficient. 

The  displacement- squared  function  * has  been  used  to  rep>resent 
the  total  damping  caused  by  struct  iiral  joints  arid  internal  materia!  JampiTig. 


It  is  represented  by  the  specific  damping  energy  whi'-h  is  the  area  within 
tile  stress-strain  hysteresis  loop  of  the  material.  Tiierefi.'re  tti''  specifi  - 
dfimping  energy  is  a function  of  the  stress  levtl  and  is  indt(  endi-.’Ht  . f t !.■ 
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shape,  stress  distribution,  and  volume  of  the  material.  The  expression 

for  the  specific  damping  energy  is  C , where  the  constant  n = 2.0  to 

3.0  for  the  usual  stress  levels  encountered  in  both  ferrous  and  nonferrous 

materials.  A rather  con^)rehensive  list  of  the  values  of  n are  given  on 

page  36-35  of  Ref.  4 . Although  the  above  equation  is  rather  general,  it 

14 

IS  impossible  to  represent  a general  hysteresis  loop. 


CHAPTER  II 


LINEARIZED  ANALYTIC  INVESTIGATIONS 
A.  Survey  of  Methods 

For  a general  damping  function,  Eq.  (l)  is  nonintegrable  and  ex- 
perimental and  analytical  methods  of  solving  the  equations  must  be  con- 
sidered. The  use  of  an  analog  con5)uter  is  not  desirable  because  of  the 
inherent  inaccuracy,  particularly  when  nonlinear  elements  are  used; 
therefore  analytic  procedures  are  considered. 

Several  analytic  procedures  of  solving  the  equation  customarily 
are  considered:  separation  of  variables,  representation  of  the  solution 

by  a power  series,  methods  of  generating  functions,  method  of  successive 
approximations,  and  representation  of  the  solution  by  finite  integrals. 

In  general,  the  separation  of  variables  is  not  possible  and  the  solution 
by  Laplace  transforms  apply  only  to  linear  differential  equations  over  a 
small  interval  of  a region.  The  other  procedures  are  generally  effected 
by  assuming  a homogeneous  and  a particular  solution.  The  homogeneous  solu- 
tion is  the  solution  of  a freely  vibrating  system  and  the  particular  solution 
is  included  to  represent  the  foi-oing  of  the  system.  As  a result  of  these 
assumptions  the  solutions  are  interpreted  as  a superposition  of  the  linear 
assumptions.  Thus  the  result  has  real  sense  only  for  linear  systems.  The 
response  for  light  or  no  damping  is  explained  as  the  superposition  of  sub 
or  super  harmonics.  The  response  for  the  damped  case  is  distorted  by  phase 
and  amplitude  variations  which  these  linear  assumptions  do  not  describe. 
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r. 


Since  the  existing  methods  are  based  on  linear  concepts,  Mii.or-ky  , 

12 

Stoker  , and  others  have  suggested  that  completely  new  concepts  might  be 
developed  to  treat  the  nonlinear  problems.  No  specific  suggestions  are 
made;  however,  Stoker  proposed  the  method  of  matching  the  response  with  a 
driving  force. 

The  general  nonlinear  response  consists  of  harmonics  of  all  orders 
and  an  excitation  force  with  a corresponding  frequency  might  excite  and 
sustain  one  of  the  harmonics.  Therefore  a plausible  physical  explanation 
for  the  effects  of  damping  should  be  independent  of  linear  ascrn;  ’ io’.is. 

There  exists  no  general  analytic  method  for  obtait.ing  an  exact 
solution  of  Eq.  (i).  Approximate  solutions  of  a knowr.  mat  he.mat  ical  form 
may  be  assumed.  Generally,  these  solutions  present  fcrm.idable  mathemati- 

i 

cal  difficulties  if  other  than  first  order  approximations  are  assumed.  The  ] 

! 

major  difficulty  is  determining  the  basic  form  and  the  required  number  of 
terms  to  include  in  the  estimate.  The  methods  vary  somewhat,  but  the  pr  - 
cedure  is  to  assume  the  solution  consists  of  linear  Independent  funct it nr . 
l.e.  , 

x(t)  = ^^(t)  + x^(t)  + D2  •••  °n  ^ ^ 

where  the  D.  are  considered  as  constants  or  functiotis  of  the  indei  r-ndent 
1 

variable,  depending  upon  the  mathematical  procedure.  In  geiicral . x (;)  Ir 

the  basic  form  of  the  solution  and  is  the  term  that  contairis  the  initial 

conditions  for  the  entire  solution.  This  solution  is  substituted  into  tlie 

differential  equation  and  evaluated  by  various  methods. 

12 

With  the  perturbation  method  , the  solution  is  approximutod  by  a 
function  of  the  independent  variable  whose  coefficient s are  i!i  p> wers  ■ f 
a small  parameter  ,|i,  of  the  nonlinear  term 

x(t)  = x^(t  ) ■*  p Xj  (t  ) • p-  X (t  ) " ■ • • • ( ) 


This  assumed  solution  is  substituted  into  the  differential  equation, 


and  coefficients  of  like  powers  of  the  small  parameter,  p,  are 
equated  to  zero.  As  the  coefficients  are  determined  they  are  substituted 
into  Eq.  (3)  to  form  the  solution.  The  term  ) is  called  the  gen- 

erating solution  and  represents  the  exact,  solution  of  the  differential 
equation  when  p,  ->  0.  The  initial  condition 


x(t  ) = X (t  ) 
' o o'  o 


(4) 


is  used  to  evaluate  the  arbitrary  constants  that  occur  when  the  solutions 
are  formed  at  each  step. 

The  frequency  of  the  oscillation  is  usually  expanded  in  a power 


series 


ID  = (JU^  M.X’  + a X 


with  the  initial  conditi  '.-.s 


x(t  ) = X (t  ) 
' o o u 


(5) 


(6) 


applied  in  the  same  manner  as  the  steady-state  amplitude  expression  above. 
This  provides  a technique  for  eiimitiating  secular  terms  which  have  tiie 
amplitude  growing  indefinitely  with  time.  Since  the  accuracy  of  the  method 
depends  on  the  number  of  correction  terms  used  and  the  magnitude  of  the 
nonlinearity,  this  procedure  would  not  provide  ati  ac(  urate  base  for  com- 
paring the  effects  of  the  dampit:g  functions.  Also,  the  perturbation  method 
is  less  useful  where  changes  in  amplitude  and  phase  occur.  This  is  pr> - 

cisely  the  major  effect  of  dampnng  on  the  response. 

J.  10 

The  method  of  variation  cf  parameters  provides  a belter  tech- 

nique for  describing  t iie  changes  iii  amplitude  and  j’hase.  This  jrocedure 
provides  for  the  cotist.ant.s  of  t ne  assumed  sol  u‘ ion  to  be  furu-t  iiiiis  ,f  the 


T 


li 

independent  variable.  The  technique  is  customarily  Illustrated  with  ti^e 
following  solution; 

X = X(t)  sin((jut  + cp(t)).  (7) 

The  derivative  of  Eq.  (7)  is 

X = X(t)  sin((jut  + cp(t))  + X(t)  cos(out  t cp(t))  [ti  + d)(t)] 

= X(t)  sin(u)t  + cp(t))  + X(t)  m cos(o)t  cp(t)) 

+ X(t)  cp(t)  cos(uut  + 9(t)).  (8) 

This  velocity  expression  must  be  changed  to  be  conq.atiblf-.-  witn 
the  original  assumed  solution,  thus  the  additional  restriction 

X(t)  sin((ut  + cp(t))  + X(t)  cp(t)  cos(uut  + cp(t))  ^ 0.  (9) 

The  acceleration  becomes 

X = ^ X(t)  uj  cos(uJt  + cp(t)) 

= X(t)  UL)  cos(u)t  + cp(t))  - X(t)  u)  sin(ujt  -t  cp(t)) 

- X(t)  cp(t)  uu  sin((jut  + cp(t)).  (lO) 

Equations  (9)  and  (lO)  are  substituted  into  the  equations  of  motioti.  The 
equations  may  be  solved  for  X and  cp. 

Because  of  the  nonlinearity  it  may  not  be  possible  to  Sidve  the 
two  equations  exactly.  In  this  case  the  values  are  customarily  t va  1 iiat  t>d 
by  the  averaging  procedure  of  Kryiuv  and  Bogoliubov.^ 

2n 

X avg  - ^ f F(x.x,t)  cos  cp  d cp  (ll) 

o 

2n 

V avg  - J F(x,.x,l)  sin  cp  d cp.  (l.-) 

o 

Thie  amplitude,  X,  and  phase,  cp,  are  treated  an. 
the  average  value. 


lb 


'.’onstaitc  when  dt  t eriiiir.  inn 


Several  methods  of  evaluating  the  residual  error  (the  aj'ithmetic 
difference  of  the  approximate  solution  from  the  exact  solutioti)  of  an 
assumed  function  are  available.  The  most  important  methods  are  the  virtual 
work,  Ritz,  Harmonic  Balance,  and  Galeriken  methods. The  approx i.mate 
solution  is  substituted  into  the  differential  equation  and  the  resulting 
error  is  defined  as  a function  of  the  independent  variable.  Since  the 
residual  is  not  a direct  indication  of  the  difference  between  the  assumed 
and  exact  solution,  the  evaluation  of  this  error  has  been  treated  by 
various  means. 

The  virtual  work  and  Ritz  methods^^  require  the  average  enej-gy  t. 
be  zero  over  the  interval  in  question  by  specifying  the  instantaneous 
virtual  work  to  be  zero.  When  the  virtual  work  expression  is  derived  from 
Lagrange's  equation. 


J e(x,t)  Z 6D  dt  = 0 , 


where  x is  the  assumed  solution,  6 is  the  first  variation  of  x with  respect 
to  D.  The  ex{  session  for  the  integral  to  be  minimized  contains  a ter.m  that 


specifies  trie  mariner  in  which  the  p.arameters  of  the  apprcxin.ate 


iTi’*  -.M>  n 


a re  t r = 


;.  The  ; ra  -ticai  value  of  knowing  the  manner  ..  f det>  r- 


1 : has  net  been  demonst. rated . 


D.  a!>'  -.lel’ficients  and  the  X.  ar-e  oi’"^hcg.  i.a.  V’-:- 
; 1 


r ' * * ■■gi*at  1 ti . 


.tin.  as  deter.mined  by  the  Ritz  [ !-■  cediure . is 
I'.-.i  fincti'TiS  that  car.  be  minim.ized  by  the  ititegrai 


I ( ' -D.  (t  )dt 


See  page  9.. . 
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In  Eq.  (2)  the  X^(t)  are  crthcgor.al  f^nctior.;  over  the  integra- 

3-' 

tion  interval.  The  D.  are  .nini-Tiized  by  requiring  - 0.  This  gives 

i 

n algebraic  equations. 

The  principle  of  harmonic  ca*ar- j-equires  the  approxi.mate 
solution  to  be  adjusted  to  satisfy  all  terms  of  the  fundamental  frequency. 
This  procedure  can  be  used  for  the  same  class  of  problems  as  the  variation 
of  parameters.  The  periodic  solution  .must  be  expressed  au  a Fourier  .jeries, 
i.e.,  the  solution  is  a periodic  function  in  time. 

In  Eq.  (2)  the  X^(t)  are  orthogonal  circular  furictions  anu  • r.e  T.. 
are  functions  of  time.  The  residual  is  formed  by  substituting  t.ne  s 
tion  into  the  integrals  of  the  Ritz  procedure 
2n 

J*  s(t)  cos  (utd(ujt)  =0  (15) 

o 

2tt 

r s(t)  cos  nuutd(uut)  = 0.  (16) 

o 

10 

The  Gaieriken  method  is  the  same  as  the  Ritz  metliod,  except, 
that  the  residual  is  minimized  by  the  principle  of  least  squares.  This 
criterion  is  more  appealing  rnathematicaJ  ly;  however,  the  nu-mber  of  al- 
gebraic equations  obtained  is  twice  the  number  ..obtained  by  the  Ritz 
method,  less  one.  The  residual  is  given  by 
t2 

J = / e^(t)  d(t).  (1-) 

The  major  difficulty  of  assuming  a functionai  f rm  f r t.he  dis- 
placement and  substituting  its  derivatives  in  the  differential  equati'.n 
is  the  error  in  the  wave  forms.  No  mathematical  pr'ced.re  l.as  beer.  f.  und 
that  will  adeciuatrly  describe  the  solution.  F.'r  •a.xs  f "j.ciu.."  da.mi.i;.g. 
excellent  results  have  been  . t ‘ ulned  for  the  lispla. . r .id  te 


noted  that  the  same  difficulties  exist  for  distributed  systems. 


application  of  Tscliebyscheff  polynominals  to  the  equation  of  motiu;.  „f  a 
fundamental  system  yield  practical  results  for  a lightly  darajied  it.fi 
medium  only.'*’^ 

At  any  rate,  the  above  methods  yield  analytical  exjjressions  that 
are  relatively  good  for  displacements,  fair  for  velocities,  and  perhaps  no 
good  at  all  for  accelerations.  Accuracy  in  the  derivatives  can  be  obtained 
only  by  an  accurate  description  of  the  displacement.  The  mathematical  dif- 
ficulties of  deriving  accurate  expressions  becomes  imjjractical  for  knovrt. 
displacement  wave  for.ms.  Since  the  nonlinear  problems  preserit  unknowt;  con- 
I ditions,  the  most  expedient  method  for  evaluating  the  general  characterist i.'S 

I 

i]  of  a system  is  the  procedure  that  would  yield  results  of  equal  accuracy  f-..r 

^ all  damping  cases. 

’ Basically,  the  numerical 'procedures  consist  of  stepwise  integra- 

, tion  of  the  equation  by  estimating  the  value  of  the  dependent  variables  at 

j each  step.  These  values  are  substituted  into  the  differential  equation  to 

I determine  the  accuracy  of  the  calculations.  The  methods  of  choosing  the 

■ values  of  the  variables  to  substitute  into  the  equation  range  from  clioosi!ig 

[ the  values  from  a random  table  to  more  efficien"  methi'ds  fiu-  analyti  -ally 

! estimating  the  solution.  The  most  efficieiit  moMK'ds  are  based  on  T'ly or 

! 

' development  about  a point;  therefor*.,  th*-  ac-uracy  is  a fur.-tur.  * f ttie 

order  of  the  development  and  the  st.ep  size.  Th.o  R'utig* -Kut  t a**  formulas,  .an’ 

a direct  extension  of  TayKu-'s  .leve Loj.ment . The  s.  . ut  i*  -a’:  i «*  ■■  !;t  inu*  d 

I U 

■ by  the  Adams- M'ul tot;  method  wl.i'.-h  weirhs  previi-js  sKpo,;  t,  i*  t*  rmlv  t!,* 
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nonlinearity  of  the  problem.  An  estimate  of  the  accuracy  of  these  methods 
can  be  made  by  comparing  the  analytical  solution  to  experimental  data  or 
to  a more  accurate  integration  procedure. 

The  inclusion  of  both  sine  and  cosine  terms  is  necessary  to  de- 
scribe the  solution  of  a damped  system.  It  generally  is  reasoned  that 
damping  causes  phase  distortion  in  all  nonlinear  problems  as  it  does  in 
the  viscous  case.  Certainly  curves  with  large  variations  in  wave  form  and 
period  can  be  described  mathematically  with  this  technique.  It  Is  impossible 
to  know  which  terms  of  the  series  to  include  in  the  estimate  of  a solution. 

In  general,  the  transient  solution  would  contain  more  terms  thati  the  steady- 
state  solution,  since  amplitude  and  phase  variations  are  greater  for  this 
condition.  This  makes  the  evaluation  of  the  phase  angle  variations  for  the 
damping  case  even  more  important. 

Since  the  Coulomb  equation  is  not  considered  strictly  nonlinear,  it 

is  appropriate  to  discuss  other  linear  procedures.  The  most  widely  used 

7 

method  of  treating  linear  systems  is  the  equivalent  energy  method.  The 
equivalent  energy  method  is  based  on  equal  energy  dissipated  in  both  linear 
and  nonlinear  dampers.  The  use  of  the  equivalent  damping  factor  in  the 
viscous  solution  gives  a direct  evaluation  ol’  the  nonharraonic  c.hararter  of 
the  motion. 

Tlie  author  feels  the  reason  that  very  little  is  kiu  wn  abouu  effec- 
tive damping  on  dynamical  systems  is  the  la-k  f kn.  wledge  f the  natural 
response  of  the  system.  Certainly  t he  { lajoi-  ova:  .a*  i*  n f a d.-imping  fun^t  ic^n 
and/or  coefficient  could  not  be  made  if  the  natura.  resj'i  i..:e  .1’  the  system 
is  unknown.  Comparison  of  peak  ami  lituder  for  t.he  meacur''me!it,  ■ f ia-mj  Inr  ' - 
efficients  is  not  sufficiently  acurati  i',  r exact  desifti  roqui  rem<  r . A 
detailed  study  • -V  the  motion  wav  f rm  and  pha.:c  variat.i.  n ov  i-  a .•  ■' 

motion  should  reyeal  s.  me  iiit'/rec  itu'  de.:i.'n  informal 


B. 


The  Eq^uivalent  Energy  Method 


The  energy  method  has  been  used  to  approximate  the  response  of 
single-degree-of-freedom  systems  subject  to  nonlinear  damping  forces. 

This  approximation  is  effected  by  replacing  the  linear  dampirig  co- 
efficient in  the  linear  solution  with  an  equivalent  nonlinear  elemen* , 
whi.h  is  determined  by  integrating  the  basic  expressioti  for  the  energy 
dissipated  per  cycle  in  the  nonlinear  damper.  This  equivalence  is  based  cn 
the  assumption  that  for  one  period  of  steady- state  mcition  the  energy  dies 
pated  in  the  iinear  damper  is  the  same  as  the  energy  dissij.'ated  in  the  not. 
linear  damper.  The  basic  assumption  is  that  the  mtjtion  is  described  by 
the  forcing  function. 

If  the  general  expression  for  the  nonlinear  damiping  force  is 
f(x,x  ),  then  the  expression  for  the  energy  per  cycie  for  tiie  solution 

X = X sin  uot  is  the  integral  of  the  product  of  the  damping  force  and 

the  velocity,  x = X uu  cos  out,  over  a cycle, 

E = ; f(x,x)  F(g  f ) dx 
cycle 

or  (l8) 

E = J f(x,x)  xdt  . 
cycle 

The  energy  dissipat.ed  per  ycle  for  any  dair.ping  I'or.-e  f(x.x  ) 

may  be  obtairied  by  subs’ ituMng  ttie  damping  for.'c  itn,.  *his  .qua’  i>  ’ . 

For  the  motion  assumed  above,  the  equation  becomes 
2tt 

E = X(ju  J ^ f(x..x)  cos  ur.it.  (lo) 

o 

Till'  ri'sull.s  .-f  sul  i t 'it  i ng  i lie  "Xpn'ssi.'i.s  !\  r l'(x.,-<)  :nt. 

Eq.  (19)  ''r.  sli.  vr.  in  Tai  le  1. 
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The  expressions  for  the  amplitudes  and  corresponding  phase  angles 


are  found  by  substituting  the  normalized  damping  coefficient  into  th»- 


viscous  equations  for  the  displacemenr  and  phase  angle,' 


2,3 


X = 


F 

k 


v'  \ ,1,  / fh  — 


(20) 


and 


tan  0 


2 r !L- 
ti;n 

j -i^y 

n 


(21) 


The  final  expressions  are  shown  in  colamns  4 and  5 of  Table  1. 

It  should  be  noted  that  the  identical  solutions  could  have  been  obtained 
with  the  Ritz  two-term  procedure. 


C.  Ritz  Method 

The  Ritz  Method^  is  a convenient  mathematical  pj’oceduj-e  for 
approximating  the  solution  to  Eq.  (l).  It  is  relatively  easy  to  include 
noniinearities  in  any  or  all  of  t.he  terms.  The  solutions  for  the  drunplng 
cases  considered  are  giver:  for  co.mparisorj  with  the  Equivalent  Ener’gy  Mi-th,  J. 

A two-term  approximation  of  the  solution 

X = X sin(ajt-0)  (—  ) 

is  assumed.  The  derivatives  of  the  solution  are  substituted  int'.  Eq.  (l4). 
These  integrals  replace  the  approju-iate  terms  of  Eq.  (l). 


TABLE 


'X  J*  sln((jut-0)  sin  (jut,d((jut) 


• I 

+ ^ r f(x,x)  sin  u)td(ojt) 
M w 


+ r X r sin((Dt-0)  sin  uutdfujt) 
M V 


- — J sin  (JUt  sin  (jutd(o)t)  = 0. 


Let  a = ujt-0, 


ujt  = a + 0, 


d(uit)  = d a 


The  first  and  third  integrals  are 


2 J*  sin(a)  sin(a  + 0)  da  . 


Using  the  trigonometric  identity 


sin(a  + 0)  = sin  a cos  0 + cosasin  0 


and  4 sin  2 a = sin  o cos  a 


the  integral  becomes 


sin^  a cos  0 t i sin  2 a sin  0 ] a a 


Integrating  and  substituting  integration  limits 


2 I — cos  0 


] = ^ 


cos  0 . 


I 
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The  integral  of  the  fourth  term  is  tt.  Thus  one  algebraic  equation  in 
X and  6 is 

2n 


-uu  Xn  cos  0 + — J*  f(x,x)  sin  ujtd(ujt) 


k F 

+ - Xrr  cos  6 = — rr  . 

m ra 


Using  the  other  orthogonal  integral,  another  equation  in  X and  6 
is  obtained: 

2tt 


-ti)  Xtt  sin  0 + - f f(x,x)  cos  uutdujt 
m J V ? 


+ - Xtt  sin  0=0. 
m 


Solving  the  above  equations  for  X and  0 yields 

' (if)' 

n M uj 


and 


tan  0 = 


ST  ''(>'•“> 

n 

i-(-  / 

n 


(23) 


(2h) 


where 


2n 


F(X,aj)  = J f(x,x)  sin  a d a 
o 

D.  Presentation  of  Anal.ybic  Results 


(25) 


The  result  of  substituting  the  expressions  for  f(x,x)  into  Eq.  (l‘,)  and 
integrating  for  F(X,0))  is  showi!  in  column  3 of  Table  1.  The  expressions 
for  the  displacement  and  phase  angle,  as  given  by  Eqs.  (23)  and  (24)  are 
shown  in  columns  4 and  5 of  Tabie  1. 

The  general  equations  for  the  maximum  steady-slate  disjU acement 
are  given  in  col  imm  2.  Table  d.  The  equations  that  are  I'btained  by 


ALGEBRAIC  STEADY- STATE  EXPRESSIONS  USED  IN  COMPUTER  PROGRAM 


SQUAPP] 


1 


2U 


substituting  the  fixed  parameters  of  the  problem  (see  Table  2)  into  the 
general  equations  are  shown  in  column  3 of  Table  2. 

The  position  of  maximum  frequency  response  as  a function  of  damp- 
ing ratio  for  the  damping  cases  can  be  determined  by  differentiating  the 
equation  of  motion  with  respect  to  the  frequency  ratio  and  equating  the 
result  to  zero.  The  value  of  the  maiximum  amplitude  is  then  determined 
by  substituting  this  frequency  ratio  into  the  amplitude  equation.  In  all 
cases  exce^’t  the  velocity- squared  case  the  separation  of  the  vai'iables 
after  the  differentiation  can  be  accomplished.  The  results  are  shown  in 
columri  H of  'latie 

The  choice  of  damping  coefficients  to  substitute  into  the  algebraic 
expressions  in  Table  2 are  not  obvious;  however,  the  equivalent  damping 
energy  method  may  be  used  to  approximate  the  coefficients.  The  expressions 
for  the  coefficients  are  given  in  Table  3-  The  coefficients  were  determined 
by  normalizing  the  energy  expression  and  substituting  the  viscous  coefficient 
into  the  resulting  equation,  the  only  assumption  being  the  equality  of  energy 
for  the  two  expressions.  The  normalized  expressions  are  giver,  in  Table  3- 

Many  authors  (such  as  Wiela^d'^^)  have  performed  experiments  by 
assuming  the  damping  coefficients  at  the  resonant  frequency  are  sufficiently 
accurate  to  determine  the  response  at  all  frequencies.  Therefcre.  the  -ia;:q - 
ing  equivalent  coefficients  are  det.errained  on  this  basis  so  tiiat  the  '-rr  r 
of  this  assumption  can  be  determined.  Also,  the  range  of  the  e^  ef  I’i c ie!:t  c 
may  be  approximated  by  this  method,  since  a value  for  the  frequei  'y  j • aK 
ampjlitudes  cannot  be  determined  beforehand.  Although  these  coeffici.c.-r  c ur- 


obviously  wrvng  for  the  exact  solution,  they  are  '.erect,  for  t.ti^'  assum-u 
linear  solution.  The  coefficients  cai.  be  evaluated  at  the  reso'i.ant  I'-aK. 


TABLE  5 


LINEARIZED  EQUIVALENT  VISCOUS  COEFFICIENTS 


TYPE  OF 
DAMPING 

LINEARIZED 
DAMPING  TERM 

EQUIVALENT  COEFFICIENT 
NORMALIZED  TO 
VISCOUS 

DAMPING  COEPTICIENT 
CORRESPONDING  TO 
ASSUMED  VISCOUS 
VALUES 

COULOMB 

^ n 

^1  u 



P^cD  nX 
^2  n 

^1  - ^ 

O.OOTT15313 
0. 3857636 
O.TT15313 
1 . 5^3063 

VISCOUS 

2 

CD  M 

R - -ii- 
^2  0) 

0.003 

0.025 

0.05 

0.1 


VELOCITY 

SQUARED 

i 

2m 

^ jto)  M 

R ^ 

5 ■ ^ 2 „ 

0)  X 

6 ^ 

^3  H CO  X 

n 1 

0.00002998179 
0.001499089  1 

0,0029^.)8it3 
0.0059<^537 

j DISPLACEMENT  , 
SQUARED 

i 

1 

2m 

, jro)  M 

i 

1 

. } 
^4=17  X 

0.231449  ! 

1 . 15725 
2.314490 
4 . 62S979 

I 
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The  degree  of  error  can  be  estimated  by  solving  the  displacement  expressions 
given  in  Table  2.  The  algebraic  equations  in  Table  2 have  only  positive 
roots  for  the  values  of  the  dissipation  coefficients  considered.  By  defini- 
tion the  magnitude  of  the  assumed  solution  is  a real  positive  t. umber. 
Therefore,  the  other  roots  of  the  equations  are  not  considered.  A computer 
program  ( see  page  123)  was  written  to  solve  the  expressions  for  the 
steady- state  displacement  and  phase  angle  equations  (column  3 of  Table  2) 
with  the  coefficients  in  Table  3-  The  values  of  the  dispiace.mtnt  and  phase 
angle  are  substituted  into  Eq.  (22),  the  assumed  solution,  which  is  sub- 
stituted into  the  differential  equation  for  the  corresponding  damping  case 
(see  Fig.  1 , p.  ?). 

The  steady- state  displacement  response  diagrams  are  shown  in  Figs. 

2 through  9-  They  are  plotted  as  a function  of  the  dajnpiiig  ratios  s.howri. 

The  root-mean- square  (RMS)  error  was  calculated  over  each  cycle  (l50  points/ 
cycle  ) . A chart  of  the  RMS  error  is  shown  in  Table  4.  The  value  of 
the  error  for  the  viscous  case  is  in  the  order  of  10  These  values  were 

checked  on  the  desk  calculator  and  found  to  be  correct.  Table  4 shows  that 
the  RMS  error  in  the  solution  for  cases  other  than  viscous  damj-ing  i-ange 
from  10'^  to  10  ^over  a cyc'e.  The  variation  for  200  points  per  eycle  gav- 
the  same  order  of  magnitude  of  the  RMS  error. 

Since  the  incremental  values  indicate  ai;  Insignificant  change  it: 
the  RMS  error,  a comparison  of  these  incremental  values  with  the  exact  s lu- 
tion  would  reveal  little  additional  information  to  add  tc  the  peak  •utijli- 
tude  analysis  above.  Also,  a det.ailed  analysis  of  itv  phase  atigle  eri\  i-  is 
not  beneficial. 


w 
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TABLE  4 

RMS  ERROR  OF  THE  ANALYTIC  SOLUTIONS  OF  THE  EQUATIONS  OF  MOTION 

The  RMS  error  of  the  analytic  solutions  of  the  equations  of  itotion 
are  given  in  the  following  p=ig.;s. 

The  equations  are  presented  in  the  same  order  as  they  appear  in 

Fig.  1.  Also,  the  same  designation  is  used.  The  dainjiiug  coefficients  are 
indicaced  in  the  second  column.  The  steady- state  am]ditude  is  showi.  in 


column  3-  The  phase  angle  (in 

radians)  is  given 

in  col  amri  4 . The 

frequer.- 

cies  are  shown  in  column  p- 

The  RMS  error  is 

shown  itj  the  last 

column.  Th 

second  line 

under  each  equation  designation  is 

the  value  of  the 

steady-stat 

displacement 

, velocity,  and 

acceleration,  respectively. 

The 

values  of  the  damping  coefficients 

are : 

Beta 

X, 

X, 

X, 

1 

2 

3 

4 

1 

0.0077153 

0 

0 

0...' 

2 

0.1 

0.005 

0. 00003 

0.  iUbol.- 
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The  iJteady-state  response  diagrams  from  the  comi'U*- r progi’a;:.  are 
given  in  Figs.  2 through  9-  damping  coefficients  used  were  caiculat'.  d ! 

by  the  equivalent  energy  methods,  assuming  the  viscous  coefficients  woluu  f 

give  a complete  response  diagram  above  critical  damjjing  and  an  output  J 

amplitude  of  unity.  These  assumptions  are  seen  to  give  excellent  cover-  ! 

age  of  the  response.  | 

The  response  curves  can  be  used  to  study  the  effects  of  da.mj)it.g.  I 

The  variation  in  the  functional  relationship  as  the  damping  cotf ficient  f 

is  increased  can  be  estimated  when  compared  with  the  nu.merical  solutf’i:  ^ 

in  Chapter  IV. 

The  values  of  the  damping  coefficients  used  for  pl-ttirig  tlu- 
following  response  diagrams  are: 


Damping 

Velocity 

Coefficient 

Coulomb 

Velocity 

Squared 

Squarod  | 

BETA 

1 

0 

0 

0 

BETA 

2 

0.10000 

0.01000 

0.00050 

0.50L>0v 

BETA 

3 

0. 30000 

0.03000 

0.00200 

1 . 90000 

BETA 

4 

0.50000 

0.05000 

0.00400 

2.50(*00 

BETA 

5 

0.70000 

0.08000 

0.00600 

BETA 

6 

0.78000 

0.  lOOOC 

0.00800 

4.500 '0. 

The  frequency  ratios 

for  the  plotted 

p lii.t.s  ap  : 

OMGP 

= .10180 

OMGP  - : . -p'.' 

OMGP 

= .PO359 

O.ttlP  - l.M  ‘ , 

OMGP 
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OMGP  - ; . 1 • 
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OMGP 

- .40719 

OMil  - i.-.^.;- 
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OMGi 

.01078 

OMGi  - . 

OMGP 

- .71G58 

0M\  II  - ^ . 

OMGP 

.Hi438 

OMIP 

OMGP 

.91*38 

OM  IP  - . . i •,/ , 

OMGP 

1.017G7 

OMil  .. 
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OMGF 

--  1.11  "77 

1 >M  II 

OMGP 

- 1 . '.‘I‘i7 

OMGI  . G.  ■ . . 

Wh>-'  P 

OMGP  - — . 

'Jj 

' r»i  1 1 . . ♦ 4 » ' 

" 

f! 

L.  
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FIGURE  4 

LINEARIZED  RESONANCE  CURVES  OF  SYSTEM  WITH  VISCOUS  DAMPING 
AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  6 

LINEARIZED  RESONANCE  CURVES  OF  SYSTEM  WITH 
LOCITY-SQUARED  DAMPING  AS  A FUNCTION  OF  DAMPING  RATIOS 


FIGURE  7 

LINEARIZED  PHASE  ANGLE  CURVES  OF  SYSTEM  WITH 
VELOCITY-SQUARED  DAMPING  AS  A FUNCTION  OF  DAMPING  RATIOS 


FIGURE  8 

LINEARIZED  RESONANCE  CURVES  OF  SYSTEM  WITH 
displacement-squared  damping  AS  A FUNCTION  OF  DAMPING  RATIOS 


FIGURE  9 

LINEARIZED  PHASE  ANGLE  CURVES  OF  SYSTEM  WITH 
DISPLACEMENT  SQUARED  DAMPING  AS  A FUNCTION  OF  DAMPING  RATIOS 


CI-IAPTER  III 


DItlEIJSIORAL  MALYSIS 

A.  Introduction 

Dimensional  analysis  of  any  physical  problem  can  be  utilized  to 
great  advantage,  the  only  requirement  being  that  one  have  sufficient 
knowledge  of  the  physical  variables  to  apply  the  procedure.  The  require 
merits  do  not  demand  a con^ilete  understanding  of  the  functional  relations 
but  rather  a knowledge  of  the  most  important  physical  variables.  Full  u 
of  dimensional  analysis  has  not  been  made  in  vibration  pi'oulems.  It  is 
believed  that  the  section  on  dimensional  analysis  is  unique  and  indi'-’ate 
the  general  conclusions  that  can  be  drawn  fi-om  atialysis  in  the  case  oi'  a 
damped  system.  To  the  knowledge  of  this  author  no  one  has  [ resentec.  an 
analysis  of  this  type. 

The  uniqueness  of  the  analysis  is  the  ap;i  1 icati^':.  of  di.monoh  :.a. 
analysis  to  the  dimensionless  ratios  to  use  in  the  resp.n.se  .iia.n-.cr..  Tn 
analysis  includes  the  definition  of  tiie  eri;i-al  ix’q  In.  rati  .’.  T.'.Ic 
allows  significant  general  concl  usiotic  in  t--  iriw.  fr  -U'  • .r 

without  knowing  the  exact  functional  i-r • 1 a t i t.;nn  i i • n-  '.  ailaih;;.  A.  c 
the  fundamental  linearization  of  the  syst’’::;  -an  !•  .ce  i • int.  rtre-  • ;.i 
results.  For  examj-'le.  a c>  ,,'npar  i son  of  t,he  ••qilva.'-n;  Iira''ir  .''e  . .a’ i onc.n  I 
and  the  dimetisionai  reiations  reveals  a •i,-f'''me>.t  ■ f ‘ h<  rec...'c. 

The  exai’t  range  >->f  this  u'reemetif  is  itUicatei  it:  t.h'  . 'ti  tiie  a.na 

of  tlie  results. 


U-i 


B. 


ResuJts  Obtained 


Considering  the  single-degree-of- freedom  system  shown  in  Fig.  1. 
the  following  physical  variables  corapLetely  describe  the  syst-nm  behavior 
(expressed  in  terms  of  the  units  force  F,  length  L,  and  time  T). 


M - 


FT“ 


F - F 


X - — 

T 


X - L 


B iS.  , • F.  • f-  . 

If  ■ L 


It  should  be  noted  tha1  all  of  the  variables  aj  i’. 

equation  of  mvition.  A dimensionally  himiogenei  u.-  • i .v  ; ■ 
variables  are  indepcndetit  aeeurding  ti.  the  bu -k  1 nmar. 

Any  relationships  that  tire  deductible  I'l-  .m  v- ■ 

consequently  dimenoiotial  )y  hemi.'genceus . Tii-'r*'  '.'•ris*.:  a 
of  the  form 

By  virtue  ef  Buck.  .•  tii-.-'.  r-'m.  <f.e  l\'dl«wi:.i' 

s in  b I'endent.  ••ath.  s ••.xl.'ts: 


' . * . 1 1 


1-1 . ■ ! .at 


of  dimensionles 


uc  X ^2 


'JO 


2^ 

F 

k 


3. 

F 


'X 


■ M 


F 

k 


9;/ 

k" 


Dimensional  ratios  other  than  the  ones  above  lo  exist;  wt 
the  ratios  shovrn  are  the  most  useful  beeause  of  their  i;;aej  • n a ;. 'v  ^ 1 
amplitude.  A compai’ison  of  the  viscous  ratios  and  i he  v.d!.'-  r i-ati.  c 
Table  5)  can  be  used  to  interpret  the  general  areas  ■•■f  th>.-  la  S]  cns<-  ui 
General  conclusions  concerning  the  various  ai-eas  ^ f ti:-.  i-'  s(  , n 


diagram  may  be  draw;,  if  it  can  be  shown  that  the  ratio's  ar>/  it.de; '‘nd'  :. 
The  dimensionless  ratios  for  C.  as  show;,  in  Table  1.  it.uicate  the 
ti-end  of  the  functicm  if  C._,  is  independent  of  the  variai  les  ■■■ 'tis i n -f' 
the  case  of  viscous  dam;!ing  (C^  is  equal  to  is  kt;  w;.  *.  -.ttue 

dent  of  the  other  physical  variables.  If  the  sam-;  inuet  it.U'.e.ei  -an  i'- 
established  for  the  other  cases  cf  dampitig.  theti  get.eral  a -t.c  I usi  Jis 
be  drawn  concerning  the  resj'e'tise  diai'rair.s. 

Upon  close  ins;iec'tii’n  oi  tli'-  rat  ios  sls.  w;  It:  Tal  i'  . ' !.■  sit!.; 
of  the  express  1 ons  f'  r C wsull  sus.gest  that  t.h<-  rat  ..-s  ;i.  ..  . . 'ir- 

critical  damping  ratios.  The  b mit.ant  ihyslca;  va.'-iat.es  1:,  •i."  ■.  1 s ■ 

critical  case  are  the  mass  atid  s;  ri’.r,  ■ '.s;  at.t  . Th"  'it.hii-  v I • ..  ••a' 

may  be  obS'-rv.  i a;:  i >miriat  it.;'  the  s;  rii:g-  • M r.  ; ; ft  i.  ■ . .■  : c\.\  . 

contr-c le,t  atid  t h.-  - '.i.t  >•■  '.b-!  p t-*  1 f • nasr-x".. 

genera!  .’.'iicl  tsiens  'at;  !•  p 'i  • ■ • ••  i-?. • • • r i'.  ' a.  - . . 

exam;le.  *!ie  S[‘ri;.g  ■ ','r  1 i t"  i '■  ''  ''■*  r-  • I sc-  ; , -i,  t'a-.  ■ ■■ 

I'jare  t-  ta;:.;  inc,  w i ; : vtir;,  ; 1 p ' : a.’  • n • i-  ' v 

x' 


at.  i X . 


( 


TABl^  5 

POSSIBLE  DIMENSIONI^SS  DAMPING  RATIOS 

I 


DAMPING 

FUNCTION 

CRITICAL 
DAMPING  RATIO 

SPRING 

CONTROLLED 

RATIO 

— ■ 

MJV3S 

GCNIR..LIED  , 
RATIO 

1 

COULOMB 

C F 

X 

C Xk 

X 

C,;-Ld'X  I 

C \ Mk 


VISCOUS 


C k 

X 


A suiniaary  of  the  tioj-n.a.  i : I'j-v;;  . ••  ■ • ffl  '1  .*  ■-  f r i ’ i-e 

equivalent  energy  method  ana  i'-:,  . . r st  1 !.■  rh  w.i  In  Tal.e  . 

The  coefficients  are  norraalizei  • ’■  '..s-  -hr.*  . -niUt  the 

equivalent  dajaping  coefficient  f i'  ■ * n -act  -fj.;.  ; ■ ax-*- i . The  tiuinex-j- 

cal  value  of  the  damping  coefficients  f r '_rc  vax-1  ■nsec  are  shewn. 
These  values  are  calculated  at  tile  I'es  •.■a:.:  j»  a.'i  f i-adlans  / 

second.  These  values  indicate  the  range  f d’axsj  ii.g  •nef I’i sients  at. 
resonance  and  are  calculated  disregai'.iir.g  any  aej e.na..!. -m  ,,n  and  acsu;:.- 
ing  an  amplitude  at  resonance  ol'  one. 

The  analytic  methods  of  investigating  the  resp  use  ef  t!i.  uatsj  ins 
cases  do  not  give  one  a clear  and  e.xplicit  feeling  f . r tile  :;.t  ins  f 
stability.  The  energy  metliod  suggests  the  correct  inteiq  retat  i- t*  {ia'^ 
on  the  meaning  of  stability.  The  procedure  is  tc  saisulato  tiie  equivai-Tit 
dampiing  coefficients  for  the  vari->us  cases  and  to  Ci’ir.pare  thesi.'  ••.■ffi- 
cients  with  the  viscous  case.  Tiiis  metliod  would  shew  the  ar’ as  _f  ‘i-- 
response  diagriuri  that  are  underdumpe  i and  'verdarr.]'eJ.  Tims,  a .'cnora. 
Interpretation  -.d’  stability  in  t.  I'mc  . f unie>-da:r.i  iiu  aV'l  i-.iamplng  is 
possible. 

The  equivalent  enei-s.v  soef t'ic ients  f t-  the  .-ases  as  wn  in 
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For  velocity- squared  damping,  the  magnitude  of  the  daraj.  ing  co- 
efficient over  the  frequency  range  considered  indicates  that  the  response 

SjjF 

is  stable  for  ail  values  of  -rf^  between  0.6  and  2.0  (predicted  by  equal 

IVIK 

energy) . 

For  displacement- squared  damping  the  coefficients  are  seen  to  vary 

considerably  with  frequency,  i.e.,  a stable  system  caiinot  be  consti'ucted 

unless  the  damping  is  varied  with  frequency. 

For  Co'uiomb  damping  the  coefficient  varies  aj q rcximately  as  f-  r 

displacement- squared  damping,  i.e..  the  response  diagram  slu  ,.ld  .-entain  a 

subresonance  above  and  below  the  natural  frequency.  Stable  resp.tise  shoula 

3, 

occur  for  values  of  between  0.39  and  1.5  (predicted  by  equal  energy). 

F 

The  principle  of  similarity  was  co.mbined  with  ttie  equivalent 
linearization  technique  to  interpret  the  response  of  the  damped  systems. 

The  dimensional  analysis  of  the  problem  revealed  some  of  the  ways  in  which 
the  constants  enter  into  the  solution.  A’thougii  mathematical  derivations 
are  frequently  attacked  with  this  general  method,  the  inter]  i’ctation  . f th' 
results  is  unique  in  that  ttiey  have  not  been  applied  t..  t.he  vit'-ati.n  ;:u;i - 
ing  cases. 

If  tiie  problem  is  solved  f.  r the  case  in  which  these  e ti.’tant  .■  ai-' 
specified  numerical  values  (i.e..  unity),  the  .u  luti  !i  may  i • im-ae  : i at  • ly 
generalized  to  cover  the  case  in  which  the  constants  havi'  ari  itrary  va.  ••  . . 
These  solutions,  along  with  the  proof  of  indepcttviciiv’c . ar'.  giv'-n  in 


Chapter  VII. 


CHAPTER  IV 
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ADAMS- MOULTON  NUMERICAL  ANALYSIS 
A.  Intruductioti 

The  method  of  obt.aining  the  soiuti'in  to  oi’diMUi-y  I'irct  ui.d  c-.  ••  r.  i 
order  differential  equations  by  ttie  fc  .ufth  oi-der  Rungt--Kutta  ana  Adair.r 

Ij. 

formulas  can  be  found  in  books  giving  a aisaussi^Vi  .1'  tuuaeri  -al  agi  r xi- 
mations.  The  Adams-Moal ton  method  for  contitmation  ^.f  the  incremoiiis 
the  function  with  the  existing  slopes  is  foutid  in  most  advanced  na’ix-i'i 
analysis  texts.  The  Adams-Moulton  pro.'cd'ore  pr^viaes  a metd;.  '!  oi  pr'  - 
dieting  the  error  after  each  integration  step.  The  procedure  is  inlierejiily 
more  accurate  than  the  Runge-Kutta  method  if  tin,  re  aia  !;e‘  abrupt  disconti- 
nuities in  the  solution.  Wlien  a discontinuity  occiu's.  the  aerivative  at 
that  point  does  not  exist.  Tlie  existence  of  all  the  derivatives  of  all 
orders  up  to  and  including  the  order  of  the  solution  must  exist;  thtrwise 
the  formulas  do  not  apply. 

The  numerical  solution  to  tiie  sec.j>;d  order  ■ rdinary  ii  lT''rt-t.t  ia ; 
ecpiations  is  accompli  sited  by  redifing  Uie  equativSis  tc’  twi>  si.mui'atc.  .0 
fii’St-order  equations  aiid  solving  these  simiU t.aneous  equati.iis  will.  Mia".- 
Mi.)Ulton  quatioriS.  The  starting,  values  1'  r the  Adam.'-M  ;1 1- n ;r 
are  obtained  with  the  Rutige-Kutta  formulas  ant  tic  itiitiai  ■ nd  1 1 . . 

The  foiirt.h  ordi  r Rungi --Kill  t a r'l’iiiiUa  mai  ),■  ! ‘ain  . 'a>"  ; 


values  i;.: 


K . • K . 

. ) . . i 
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There  are  two  Adams -Mo  alt on  formal  as.  One  is  the  prediction  equation 

1 


y 


M ' h " it  [»  h ■ » F,  * 37  Fj  - 9 (.7) 


and  the  correction  equation 

h 


»1,1  ° »1  * 2t 


[9  ha  * 19  F.  - 5 F,^_J  « F 


3.  J • 


(25) 


The  addition  of  the  step  values  of  hAy^  to  the  correct  value 

of  is  carried  out  in  double  precision,  i.e..  72  bit  w^ rd  length  (all 

the  precision  available  on  the  CDC  3^00  computer).  Ail  . tlu-r  i,!9-rali.  n^j 

are  done  in  single  precision,  i.e.,  2h  bit  word  length.  This  canservis 

r "’h 

the  accuracy  of  the  above  methods  which  have  ei-r. cs  less  thaii  1 i.J' . 

(Note:  For  the  equation  used,  the  program  will  not  run  with  any  repeat- 

able accuracy  without  the  double  precision  arithmeti  ' . ) 

The  error  analysis  consists  of  checking  the  mnative  ma.qiitude 
of  the  corrected  value  at  each  step  by  tlie  following  f 

, ,,  r c F 1 ' 9 
-‘9Ln.-i  - ha  j 


error  = max 


270 


r c 


I 

j 


( ') 


where  C and  P refer  to  corrected  and  predicted,  rcrj  • ei i v : y . Equal!  •. 
(29)  represents  the  error  due  to  the  truncation  of  Eqo.  (,  ’)  an  1 (.  " ) . 

This  procedure  ['rovides  some  estimate  of  Uu  numt  er  ;■  si-'nif.- 
cant  digits  obtained  in  the  corrected  solut  ion  as  t he  '•alcu.  a'  i 'lu'  a:--' 
performed.  This  procedure  will  not  guarantee  the  a’^.-uracy  f t!;>  s.  ; 
tion  since  accumulative  err  rs  eannot  *.17  del.u't’'i.  Als.  . the  s!'"ifi  -a- 
t ions  of  the  limit  val  les  of  the  err.  r ar>  'fit  lea:.  a!i  i will  vary  with 
the  rate  of  change  of  the  lerivativ  Th"  •■xact.  s.  tail-  j.  , ;■  E).  (:) 
is  obtaine.l  if  the  foil,  wing  ii'i.tity  i.’  satisfi'  i: 

■*  f f(x.x)  * kx  - F si!i  X'  . 


Therefore,  the  determination  of  the  residual  in  the  above  equation  will 

indicate  the  accuracy  of  the  solution.  The  root-mean- square  (RMS)  value 

of  the  residua]  is  determined  over  the  steady-state  period.  This  residual 

was  found  to  be  small  foi-  a minimum  of  150  points  per  cycle.  The  minimu-m 

and  maximum  number  of  points  per  cycle  is  fixed  for  each  frequency  so 

that  the  accuracy  and  storage  space  requirements  are  preserved.  A table 

of  the  integration  steps  versus  frequency  was  compiled  for  this  puri'ose. 

-8 

The  numerical  value  of  the  RMS  residual  is  less  than  10  for  all  ttie 
equations. 

The  steady-state  condition  exists  if  the  wave  form  of  the  solu- 
tion and  all  of  the  derivatives  vary  in  a continuing  periodic  .manner, 
i.e.,  f(x)  = f(x+r).  The  steady-state  check  is  performed  once  pei'  cycle 
when  the  displacement  changes  sign  I'rom  negative  to  positive.  Wnen  this 
amjiiitude  sign  change  occurs  all  of  the  calculated  values  over  a cycle 
are  checked  for  muxlimun  value  by  substitution  into  subroutine  AMV'.  Tile 
.maximum  value  is  ctox-ed  and  compared  to  the  previous  maximum  values,  as 
showT.  in  Fig.  10.  Tlie  amplitudes  are  compared  with  the  I'olUwing  formula; 
Dx ( 1 ) - Dx ( i - ! ) . 0.005 

If  the  amplitude  check  falls  within  the  defined  limits  tlie  peril- dicjty 

is  checked  by  the  following  formula: 

lb  (i ) - hs(l-I ) - 0.005 

Bz(i) 

If  the  motion  does  not,  .;r,'ss  the  t.lm.o  axis,  a si.  r-it'e  limit  st  i .•  th<' 


progi'um  and  the  data  is  j rinti- i and  j lotted.  (’•'-■'i  B*  th  ji  si*  iv<-  anu 
negativ"  values  uf  t.he  ■emplitaf  and  p"!'i"d  ap-  t ■ ■ rnii  t t • ■ i . ) 


5^ 


Wlrien  the  continuing  periodicity  of  the  wave  form  is  established, 
the  steady-state  and  transient  data  are  printed.  All  of  the  calculations 
are  self-explanatory  except  the  energy  dissipated  per  cycle  and  the  phase 
angle. 

The  actual  damping  energy  per  cycle  is  calculated  with  the  for.m'ulas 
given  in  Chapter  III.  The  correct  expression  for  each  damping  case  must 
be  used;  therefore,  a program  was  written  for  each  danipiing  function. 

The  areas  under  the  cuc’ves  are  determined  by  Simpson's  l/'l  rule. 

The  justification  for  using  the  1/3  rule,  rather  than  a higher  approxisia- 
tion,  is  the  inherent  accuracy  of  the  calculated  values  (see  section 
above).  Also,  the  number  of  calculated  points  per  cycle  is  Uut  neces- 
sarily an  even  number  as  required  by  other  procedures  such  as  the  3/8 
rule.  For  an  R14S  error  of  10  the  value  of  the  energy  caiculatio.n  is 
accurate  to  four  places  (a  relative  error  of  0.01'^  for  tlie  viscous  case). 

The  phase  angle  between  the  forcing  function  and  the  dicplacemeTit 
amplitude  is  calculated  at  the  zero  crossing  and  at  the  peak  a.mplitudes. 
This  is  done  by  comparing  the  angular  travel  ^f  tlie  forcing  function  to 
the  resulting  displacement  travel  with  time  equal  zero  as  a rt-fi.Te:;cc . 
The  phase  angle  at  the  zero  crossing  was  used  as  ar»  aid  in  reducing  ih'- 
data.  Unless  othei-wise  specified,  all  refei'oi.ce  to  phaso  ai.cie  will  an 
phase  angle  at  the  positive  p)eakc. 

The  values  printed  under  transient  respv-nse  arc  cy-lio  va: 
and  are  printed  without  regard  to  sign. 

A .■.,;nplete  time-hisioiy  , q'  the  1 i cpl  a ■en,c;U  . vj . '.-i  1 . ati  : a.--,  hi-- 
atlori  is  plotted  at.  t.h>'  eti.i  of  i aei.  rui;.  Ev-  r,v  t'-ntii  valu'-  1’  t h-  ia*  a !. 
plotted  I't.  a »-in.  ny  "-i!..  axi.:  'In.  seal--  ■ !'  whi-li  is  i-  t'.'rndin  u t\\- 


subroutine  Largest.  The  original  program  contained  a routine  to  plot  the 
complete  steady- state  and  transient  response  diagrams;  however,  this  was 
found  to  be  i.mpracticai  because  of  the  computer  time  required  (eight  hours 
for  viscous  damping).  Ttie  program  is  given  in  the  A]'pendix,  j . 128. 

B.  Numerical  Results 

The  computer  program  described  above  is  used  to  solve  the  equation 
of  motion.  With  the  exception  of  additional  values  of  the  dampdng  co- 
efficients, the  parameters  (shown  in  Table  3)  are  used  for  all  damjdng  case 
The  equations  for  Coulomb,  velocity-squared,  and  displacement- squared  are 
presented  since  the  viscous  ease  is  given  in  Chapter  II. 

The  viscous  case  was  used  as  a base  for  comparing  the  ae’curacy  >-f 

the  other  equations.  A direct  comparison  of  equivalent  energy  metliod.  the 

numerical  method,  and  the  results  from  a desk  calculator  revealea  :io  ilo- 

tinguishable  error  in  the  maximum  values  of  the  displ acement . veiooity.  a.u  : 

acceleration.  Several  different  frequencies  wiu'e  checked  on  cither  side  of 

resonance.  Also,  the  values  of  the  RMS  err.  r ure  oV  the  samt'  .magnit s..ie 
1 0 

(10  ' ) for  each  integration  step  as  it  is  ever  a cyole  .1’  mciio!!.  It  sh.  u 
be  noted  that  the  CDC  3200  co;ii])uter  has  e?jly  10  digits  of  a..'ouracy.  Tliei-.  - 
fore,  RMS  errors  in  the  same  order  of  matduluue  are  ectisidei-vd  s.iff : -i..  nt.  ;,v 
accurate. 

The  steady-state  response  eiiai'a.d.er-ist  Ics  el’  C'oul  rii.  veI.o-i:,v- 
squared,  and  displaceim?nt-square.i  damping  are  given  in  Fi  11  tiire  .gn 
The  RMS  errei-  is  in  lh^‘  ..rder  ..f  for  all  .'aovs.  TIk'.v  ■i!-..'  pA 

as  a futi'-'lLiii  of  tin.'  damping  rat  i^  s rivtui  i:i  Ci.apt.;-  III.  Tlie  ; 1 rp : a •em';.t 
phaS'.'  aii.gle,  vole, 'it;/,  a - ■odi'raf.ien.  a;.  : Ur-  ,ia;::|irig  >-g;,-  ;-  v -v  • a-- 


given.  Tlie  valui'S  [!■  ' 1 are  f.  r the 
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FIGURE  11 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY, 
AND  acceleration  FOR  NO  DAMPING, 

OJ  UJn  0.255,  U)  (jJr.  = 0.509,  AND  (M  'aJp  = 0.764 


FIGURE  )2 

time-histories  of  displacement  velocity, 

AND  ACCELERATION  FOR  NO  DAMPING, 


FIGURE  13 

time-histories  of  displacement,  velocity, 

AND  ACCELERATION  FOR  NO  DAMPING, 

CU  '(Vr  - 1.272,  co'cur^-  1.527,  AND  UJ'UJr^'  2.545 
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FIGURE  14 

time-histories  of  displacement,  velocity, 

AND  ACCELERATION  FOR  COULOMB  DAMPING  WITH/3,  F = 0.77, 
ou  Ldn  = 0-255,  cj  LUn  ^ 0.509,  AND  u 'uJn  ^ 0.764 


FIGURE  15 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY, 

AND  ACCELERATION  FOR  COULOMB  DAMPING  WITH/Q,  F = 0.77, 
OJ  'cu^  = 0.916,  1.018,  ANDcu  = 1.272 


FIGURE  16 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY, 

AND  ACCELERATION  FOR  COULOMB  DAMPING  WITH/3,  F -■  0.77, 
w = 1.527,oj  --  2.036,  ANDu;gj^  ^ 2.545 


FIGURE  17 

COULOMB  displacement  CURVES  BY  NUMERICAL  PROCEDURE 
AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  18 

COULOMB  PHASE  ANGLE  CURVES  BY  NUMERICAL  PROCEDURE 
AS  A FUNCTION  OF  DAMPING  RATIOS 


/3,/F  = 0.0077 
/3,  'F  - 0.385 
/3,/F  = 0.5 


FIGURE  19 

COULOMB  ENERGY  CURVES  BY  NUMERICAL  PROCEDURE, 
AS  A FUNCTION  OF  DAMPING  RATIOS 


FIGURE  20 

COULOMB  VELOCITY  CURVES  BY  NUMERICAL  PROCEDURE, 
AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  22 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  VELOCITY  - SQUARED  DAMPING  WITH  /33F  Mk  ^ 3.09, 
cutUn  - 0.255,  cu'UJn  = 0.509.  ANDcu'CLin:  0.764 
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FIGURE  23 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  VELOCITY -SQUARED  DAMPING  WITH/SjF  Mk  ; 3.09, 

(jj/UJ^  - 0.916,  ojtd,,  = 1.018,  ANDcu  'UJr,  = 1.120 


FIGURE  24 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  VELOCITY  - SQUARED  DAMPING  WITH  /SjF  Mk  = 3.09, 

(jj  'ojr,  - 1.272,  oj  'cUn  ' 1-527,  AND  Cl)  (jUp  - 2.545 
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FIGURE  26 

VELOCITY-SQUARED  PHASE  ANGLE  CURVES  BY  NUMERICAL  PROCEDURE, 
AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  28 

VELOCITY-SQUARED  VELOCITY  CURVES  BY  NUMERICAL  PROCEDURE, 
AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  79 

VELOCITY-SQUARED  ACCELERATION  CURVES  BY  NUMERICAL  PROCEDURE. 
AS  A FUNCTION  OF  DAMPING  RATIOS 


FIGURE  30 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  DISPLACEMENT- SQUARED  DAMPING  WITH/3,F  0.308, 
uj'ujr,-  0.25S,  cutUp  = 0.509,  ANDojLjp-  0.764 
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FIGURE  31 

TIME -HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  DISPLACEMENT-SQUARED  DAMPING  WITH /3,i  F/k^=  0.308, 
oi  LUn'  0.916,  Cl) 1.018,  AND  ajLL)^,  = 1.120 
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FIGURE  34 

TIME -HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  DISPLACEMENT- SQUARED  DAMPING  WITH /S^F  V r 1.5, 

1.018,  ojtUp=  1.120,  ANDco  'co,,--  1.272 
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FIGURE  35 

TIME-HISTORIES  OF  DISPLACEMENT,  VELOCITY,  AND  ACCELERATION 
FOR  DISPLACEMENT-SQUARED  DAMPING  WITH /34F  V 0 .5, 
cu 'ojp  - 1 .527,  cutOn  - 2.036,  ANDcu  xun'  2.545 
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FIGURE  36 

DISPLACEMENT-SQUARED  DISPLACEMENT  CURVES 
BY  NUMERICAL  PROCEDURE, AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  37 

DISPLACEMENT-SQUARED  PHASE  ANGLE  CURVES 
BY  NUMERICAL  PROCEDURE, AS  A FUNCTION  OF  DAMPING  RATIOS 
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FIGURE  38 

DISPLACEMEMT-SQUARED  ENERGY  CURVES 
BY  NUMERICAL  PROCEPURE,AS  A FUNCTION  OF  DAMPING  RATIOS 


FIGURE  40 

DISPLACEMENT-SQUARED  ACCELERATION  CURVES 
BY  NUMERICAL  PROCEDURE  AS  A FUNCTION  OF  DAMPING  RATIOS 
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the  Introduction.  The  curves  are  extended  by  interpreting  the  time- 
histories  where  the  steady- state  conditions  are  not  well  defined.  Where 
insufficient  repeatability  of  the  motion  exists,  no  points  are  plotted. 

The  curve  extensions  are  intended  to  aid  in  the  interpretation  of  the 
response. 

Typical  time-histories  for  each  damping  case  precedes  the  response 
curves.  These  curves  show  the  steady-state  and  transient  response  for 
average  values  of  the  damping  coefficients.  In  general,  the  displacement, 
velocity,  and  acceleration  curves  differ  signif icantly  for  each  case  and 
each  frequency.  In  all  cases  the  higher  derivatives  cannot  be  related 
with  the  simple  algebraic  or  trigonometric  equations  indicated  by  Eq.  (2). 
The  corresponding  time-histories  for  no  damping  are  given  in  Figs.  11 
through  13.  These  curves  are  presented  for  comparison  with  the  damped 
cases. 

Typical  Coulomb  damped  time-histories  are  given  in  Figs.  14 

through  l6  . The  displacement  wave  forms  seem  to  vary  as  expected  with 

the  stick-slip  motion  occurring  at  — =0.5  and  smooth  uniform  motion  at 

'^n 

the  higher  frequencies.  The  velocity  and  acceleration  curves  are  not 

uniform.  The  effect  of  the  discontinuities  may  occur  in  the  velocity 

= 0.255)  or  the  acceleratioti  - !.27^  and  higher  ratios).  The 

presence  of  low  frequencies  in  the  outjiut  is  most  evident  in  the  acceiera- 

tion  for  frequency  ratios  from  — = 0.764  to  — = l.Olo.  Tiie  envelope 

appears  to  be  described  by  a straiglit  line  with  sharj’  discontinuities. 

This  is  not  the  case  foi'  viscous  damping  wliere  the  envoiope  is  sinusoidal. 

The  wave  shape  of  the  Uiird  derivative  could  be  estimated  from  the  curve.'; 

uu 

of  velocity  and  accc  leral  ion.  For  r*xample.  the  velocity  for  — ~ 0.t>09 
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i 
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is  similar  to  the  acceleration  at  — =0.764,  which  suggests  the  "jerk" 

“^n 

would  be  similar  to  the  acceleration  at  — = O.SOQ. 

uu 

n 

The  velocity- squared  time-histories  (Figs.  22  through  24  ) are 
more  uniform  than  Coulomb.  The  wave  forms  are  practically  identical 
for  all  coefficients  and  frequencies  considered.  The  velocity  peaks  are 
smoothly  rounded  and  the  acceleration  peaks  are  sharp.  The  accelera+.ion 
wave  form  at  — = 0.255  is  similar  to  the  displacement  wave  form  (the 

negative  peak  being  disregarded)  for  — = 0.255  no  damping.  There- 

fore,  the  higher  derivatives  are  suggested. 

The  displacement- squared  time-histories  curves  are  given  in 
Figs.  30  through  35-  In  general,  the  response  is  similar  to  the 
Coulomb  curves.  The  motion  appears  to  have  the  stick- slip  type  motion 
for  small  dan^jing  and  low  frequencies  (Figs.  36  and  4o).  At  higher  fre- 
quency ratios  the  system  appears  to  be  undamped.  Figs.  11,  32,  and  35 

for  — = 2.545  are  identical.  Tiie  low  frequency  respionse  for  high  damp- 

% 

ing  coefficient  (Fig.  33  ) seems  to  be  dissimilar  to  the  undamjied  resp'onse. 

The  phase  angle  response  diagram  for  Coulomb  damping  gives  a 
better  understanding  of  the  motion  characteristics  than  the  di  sjii  acement, , 
velocity,  and  acceleration  response  curves.  The  pliase  angle  curves  (Fig. 
18)  show  that  the  response  has  a subresoiiant  peak  at  a froqueTicy  rati,.)  of 
0.5;  however.  Fig.  l4  siiows  that  stick-slip  motion  occurs  at  this  fre- 
quency and  damping  ratio.  The  [jhase  angle  in  this  region  undergoes  a 
smooth  change  in  slope  for  values  of  the  damj  itig  rat  i.  i-  tw'  ,■>'  0.  <8  a;.u 
0.11.  For  values  in  the  order  of  0.0077  the  syc'^’m  is  ur.stnid.  f,  r ff  - 
quency  ratios  of  0.5  throug.h  1.0.  Tlie  revrsal  iti  phas-  a*  - l.o 

i 1 i t y 


appears  to  be  instar'' aneous. 


For  frequencies  near  r, •sotniic. ■ 


inst  al 
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occurs  for  all  values  of  the  phase  angle  that  are  greater  that:  l80 

degrees.  Also,  the  motion  is  cort^tletely  stable  when  the  phase  angle 

drops  below  l80  degrees  for  the  frequency  ratios  considered.  However, 

as  the  phase  angle  approaches  90  degrees  for  the  higher  frequency 

ratios,  the  motion  becomes  unstable  again.  The  other  Coulomb  response 

diagrams  indicate  no  unusual  characteristics. 

The  response  curves  for  velocity- squared  damping  indicate  no 

characteristics  similar  to  the  Coulomb  case.  The  phase  angle  curves 

(Fig.  28)  are  smooth,  continuous  curves.  The  phase  shift  is  practically 

constant  for  values  of  the  damping  ratio  between  1.158  and  2.12. 

In  general,  the  displacement -squared  response  is  similar  to 

the  Coulomb  curves.  The  phase  angle  diagram  (Fig.  37  ) indicates  a phas“ 

shift  at  — =0.5,  and  the  same  instability  comments  apply.  The  motion 

n 

is  unstable  for  frequency  ratios  between  O.5  and  1.0  if  the  curve  drops 
below  zero  degrees.  For  frequency  ratios  above  1.0  and  damping  coeffi- 
cients greater  than  0.5j  the  motion  is  stable.  Da.mping  coefficients  in 

the  range  of  0.1  have  a completely  different  shape  that  would  seem  to 

uu 

Indicate  another  resonance  point  at  — = 2.0.  Ajain.  all  vaiues  of  the 

phase  angle  larger  than  18O  degrees  and  less  ’ha'  'K'  .tegr'-:.s  give  un- 
stable motion. 

In  general,  the  Coulomb  damj  i!.g  has  th<-  ’.a;-r  w-s’  sot.aru  jeak. 
and  the  displacement  - squared  * he  widest.  The  v<';>' ■ , -s  j .ic’ .1  resp:  t.se 
is  smoother  and  has  the  high-’St  . r.t  rg;.  1 i ss  i j a’ i ■ ■ . 

T>ie  position  of  maxim'.jn  fr<  1 i- r .'j  cs  a.-  a fs. 'tiv  '..  v 1’ 


damping  ratio  for  the  < 

damping  asa  .•  -an  1 ■ !•  ' 

' 1 t y 1 : ! 

t lit  i 

the  equation  of  mot  iori 

wi ' h t • S[  

j ■ . r’l*  . '1*  1 
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the  result  to  zero.  The  displacement,  velocity,  and  acceleration 
viscous  damped  natural  frequencies  are: 

i 

2 ^ 

Displacement  resonant  frequency  u;  (l-2Q  ) 

Velocity  resonant  frequency 

Acceleration  resonant  frequency  ^ . 

This  frequency  ratio  can  be  substituted  into  the  equation  of  motion  to 
determine  the  value  of  the  maximum  response  amplitude.  Since  the  functional 
relationship  is  unknown,  these  values  may  be  obtained  from  the  response 
diagrams.  The  results  for  displacement  are  shown  in  Figs.  4i  and  k2  f>>r 
the  viscous  and  velocity- squared  cases.  The  Coulomb  and  displacement- 
squared  cases  are  not  plotted.  The  Coulomb  response  curve  shows  resonance 
at  a frequency  ratio  of  one  and  the  displacement- squared  case  is  fre- 


quency-dependent . 
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FIGURE  42 

RATIO  OF  AMPLITUDE  OF  MAXIMUM  FREQUENCY  RESPONSE  TO  INPUT  AMPLITUDE  AS 
A FUNCTION  OF  DAMPING  RATIO  FOR  VISCOUS  AND  VELOCITY-SQUARED  DAMPING 
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CHAPTER  V 


DISCUSSION  OF  RESULTS 

A.  Discussion  of  Root -Me an- Square  Error 

The  accuracy  of  the  solutions  to  the  differential  equations  is 
based  on  the  evaluation  of  the  magnitude  of  the  residual.  The  residua^ 
is  defined  as  the  arithmetic  difference  between  the  approximate  and  true 
solution.  It  is  obtained  by  substituting  the  estimated  values  of  the 
dependent  variables  (for  corresponding  values  of  the  independent  variable) 
into  the  differential  equation;  therefore,  it  represents  an  aggregate  oi' 
the  deviation  of  the  dependent  variables.  Also,  tlie  residual  represents 
the  sum  of  both  positive  and  negative  terms  which  are  either  very  large 
or  small.  The  sign  of  the  terms  are  not  random  since  the  values  of  tlie 
dependent  variables  are  in  the  neighborhood  of  the  correct  value. 

The  deviation  of  one  of  the  independent  variables  from  tiie  exav't 
value  will  vary  from  point  to  point  over  the  cycle.  The  magnitude  of  the 
residual  will  depend  upon  the  integration  method  and  the  step- size.  If 
the  analytic  procedure  is  legitimate  and  the  integration  interval  is  in- 
finitesimally small,  the  approximate  value  will  yield  a small  residua!. 
Therefore,  analytic  methods  of  evaluating  the  error  may  be  applied. 

In  general,  the  residuals  will  be  of  either  algebraic  sign  with 
respect  to  the  true  values.  Therefore,  the  square  of  the  residual  can 
be  used  to  represent  t.he  error.  By  the  juincijile  of  least  squares,  tin 
expression  for  t.he  error  over  a cycle  is 
.f  e‘^(t)dt. 


I 
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The  value  of  the  average  error  is 

e 

rms 

or  the  root-mean- square  error.  By  the  same  reasoning,  nhe  average  error 
of  all  the  dependent  variables  is 


The  root-mean- square  error  was  calculated  at  each  integration 
interval  (not  less  than  I50  points)  and  the  magnitude  of  the  error  did 
not  change.  The  total  root-mean- square  error  over  the  cycle  was  equal  to 
the  point-by-point  value.  Therefore,  neither  the  calculation  procediure 
nor  the  equipment  caused  any  propagation  errors.  The  RMS  error  for 
viscous  damping  is  less  than  10  Also,  the  errors  in  the  steady- 

state  amplitudes,  displacements,  velocities,  and  accelerations  for  the 
viscous  case  is  less  than  0.06%.  Therefore,  it  is  concluded  that  an  RMS 
error  of  10  is  a reliable  measure  of  the  accuracy  of  the  solution. 

The  RMS  error  for  the  nonlinear  damped  solutions  by  the  linearized  pro- 
cedure  is  in  the  order  of  10  . In  general,  results  with  residuals  of 

this  order  of  magnitude  are  completely  unreliable.  It  is  impossible  to 
ascertain  whether  the  error  of  the  linear  procedures  occiu-s  in  the  dis- 
placement, velocity,  or  acceleration.  The  steady- state  velocity- squared 
displacement  and  phase  angle  response  diagrams  for  the  linearized  solu- 
tions are  shown  in  Figs.  46  and  47»  respectively.  Tlie  errors  in  the 
diagrams  are  given  as  a percent  of  the  n'jmerical  values.  Errors  (by  the 
linearized  procedure)  in  the  steady- stat.e  amplitudes  for  velocity-squared 
damping  are  as  high  as  0 per  cent.  (Note:  The  errors  in  the  linearized  prc'- 
cedure  for  the  Coulomb  curves  are  as  large  as  ‘jOO  per  cent.)  The  phase  at.gle 
response  diagrams  show  similar  errors.  A eomparison  of  steady-state 


9h 

displacement  errors  with  the  corresponding  velocity  and  accelera’,  ion 
errors  is  given  in  Table  6.  The  errors  in  the  displacement,  velocity, 
and  acceleration  seem  to  vary  in  a random  manner  with  no  particular 
pattern.  Although  these  comments  refer  to  maxim'jm  steady-state  values, 
the  instantaneous  values  vary  in  the  same  way.  The  instantaneous  root- 
mean-  square  errors  are  of  the  same  order  of  magnitude  as  the  RMS  total 
error.  The  wave  forms  for  velocity- squared  damping  differs  somewhat 
from  sine  waves,  and  the  Coulomb  wave  forms  bear  little  resemblance  to  a 
sine  wave.  Perhaps  the  explanation  for  equal  RMS  values  is  that  the 
Coulomb  case  has  many  solutions  due  to  the  variable  initial  conditions 
at  discontinuities  in  the  damping  function. 

These  comments  woald  seem  to  indicate  that  n ■ matT.er  what  liiiear 
analytic  form  of  the  solution  is  assumed,  t.he  error  would  be  large.  The 
method  of  optimizing  the  coef ficietits  and  functions  of  Eq.  (.  ) is  not  as 
important  as  having  the  correct  furctior.s. 

B.  Equivalent  Energy  and  Ritz  Methods 

The  accuracy  of  the  equivalent  energy  concept  is  investigated  by 
comparing  the  results  in  Table  4 with  the  accurate  value  determincu  by 
the  computer  program.  The  calculation  procedure  and  acciu'acy  used  it: 
the  computer  program  is  discussed  in  Section  A.  The  procedure  was  to 
choose  coefficients  at  the  resonant  frequency  based  on  equal  energy  at:J 
to  compare  the  results  with  the  values  obtained  witii  the  tiumerical  pr*  - 
cedure.  The  equivalence  of  ene^-gy  requires  the  assiun]>tion  of  the  steady- 
state  amj)litude.  If  t.he  wave  forms  of  displacement  atid  velocity  arc  sii.c 
waves,  the  coefficient  will  give  the  same  aniplit.ude  for  tiie  assumed  d;ur.i - 
ing  function. 
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Table  6 

Velocity  and  Acceleration  Errors  Corresponding  to  Displacement  Errors 

in  Fig.  46,  j). 


Damping  Frequency  Rati  .■  Velocity^  Acceleration* 

Coefficient  -d  Error  Error 


0.0014991 

1 

0.509 

+1.14 

-24.2 

1.081 

+0.89 

-13.4 

1.527 

+0.61 

-3.44 

0.0029982 

0.509 

+3.22 

-52.5 

1 . 081 

+1.61 

-22.2 

1.527 

+1.22 

- 7.7 

2.036 

+0.25 

- 2.88 

0.0043163 

0.509 

+4.6 

-75 

1 . 081 

+2.1 

-26.6 

1.527 

+1.59 

-11.3 

1 

2.036 

+0.62 

- 5.12 

0.0059964 

0.509 

+5.7 

-102 

1.081 

^2.5 

-32 

1.527 

+1.9 

-15 

2.036 

+0.95 

-7 

* 

The  errors  are  given  as  a per  cent  of  the  numerical  value,  i.e.. 

Numerical  Value  - Linearized  Valu', 

Numerical  Val  ue 


L J 


I 

I 
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The  accuracy  of  these  assumptions  for  velocity- squared  damping 
and  viscous  damping  is  shown  in  Figs.  43,  44,  and  45.  An  amplitude  of 
unity  was  assumed,  and  Fig.  43  indicates  the  assumption  is  correct  for 
the  displacement;  however,  the  velocity  and  acceleration  values  differ 

i 

significantly.  Also,  it  should  be  noted  that  the  "fit"  of  the  curves  is  ) 

J 

1 

for  an  insignificant  frequency  range.  j 

i 

Since  the  wave  foim  is  the  most  important  consideration  for 
equality  of  motion,  the  equality  of  energy  is  a good  appi'oximat  ion  for 

first  order  accuracy.  If  the  correct  amplitude  could  be  determined,  the  i 

energy  procedure  will  give  accurate  results  for  a fixed  frequency.  The  ; 

energies  are  equal  for  the  same  damping  function  if  the  wave  form  is  not 

changed.  The  only  requirement  is  the  equality  of  the  dimensional  ratios. 

The  representative  errors  in  the  linearized  steady- state  velocity- 

squared  response  diagrams  are  showi;  in  Figs.  46  and  47.  The  errors  are 

expressed  as  a percentage  of  tlie  numerical  value,  i.e.. 

Numerical  Value  - Linearized  Value 
Numerical  Value 

The  errors  are  representat  iv^e  of  the  variation  from  the  assumed 
solution.  In  those  areas  of  the  resj  ^.nse  diagrams  where  the  error- 
approaches  zero,  the  functional  variations  are  negligible.  Tire  most  signi- 
ficant deviations  are  indicated  tiy  errors  as  large  as  0 per  cent  and  as  small 
as  0.01  per  cent,  depending  upon  frequency.  T;, -pi cal  velocity  and  acceleration 
errors  for  the  displacement,  are  given  in  Table  i . The  velocity  err-or-s 
are  positive  and  the  acceleration  ei-rorc  are  negative  for-  all  fi-eiiuency 


ratios.  The  larger  errors  oeeur  at  the  lower-  frequencies  in  ail  cases. 
Because  of  the  wide  variati.m  i-f  th'-  full  •-.cverage  of  the  diagi-am;- 

was  considered  unnecessar-y  si:  a '-iU'-ite  ■■.•-ve  - ar-c  given  in  Cliapt--!-  IV. 


FIGURE  43 

COMPARISON  OF  DISPLACEMENTS  OF  VELOCITY-SQUARED  AND  VISCOUS 
DAMPING  FOR  EQUAL  ENERGY  DISSIPATED  AT  RESONANCE 
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! 
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FIGURE  46 

REPRESENTATIVE  ERRORS  FOR  RESONANCE  CURVES  OF  LINEARIZED  METHOD 
-VITH  VELOCITY-SQUARED  DAMPING.  FOR  VARIOUS  VALUES  OF  THE  DAMPING  RATIOS 
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C.  Diif.L-tisior.ai  Afialyji^ 

Demonstration  of  t.^ie  Itidependence  of  t.he  Dlmensicriless  Ka’.  Iqj 

The  dimensionless  ratios  givetj  in  Table  5 he  shovm  to  b>. 
linearly  independent  by  matrix  theory;  tiowever,  tiie  inaepeiidet.r-  is 
demonstrated  to  guarantee  that  no  mathematical  errors  have  occurred. 

The  independence  of  ttie  ratios  is  demonstrated  oy  solving  the  differet  - 
tlal  equations  for  various  values  of  the  parameters.  By  changing  one 
parameter  at  a time  the  prediction  of  the  response  is  confirmed. 

The  effect  of  changing  the  damping  coefficient  is  evid',-’.’  i,v 
studying  the  response  diagrams.  As  witii  the  viscous  case,  remai:..' 
unchanged,  and  the  damped  natural  frequencies  vary  somewhat,  depending 
upon  the  type  of  damping.  Figure  Ui  shows  tlie  v^ariation  of  the  velc.'ity- 
squared  and  viscous -damjied  natural  I'requefi'-’ies . Tlie  displacemc;.’ -rqifartu 
and  Coulomb- damped  maximum  frequency  is  giveri  in  Table  e.  In  all  cases 
the  damping  ratio  varies  directly  as  the  cliaT.ge  in  the  dairiping  CL-effi-- i-.f'.t . 

The  effect  of  changing  the  for  ifig  anii  lit..de.  mars,  and  sj  ring 
constant  is  demonstrated  by  hojding  aii  of  the  variables  ccnstaf:*  ex -i  i * 
the  one  being  considered.  Tiie  if.aeper.der.cc  f .a’!:  cf  t.he  varia't  les  i r 
verified  in  this  mafuier.  T.hi-  • fJ'c -t  -f  ciiangitig  the  above  variat  l--s  It  r 
viscous  damping  is  well  kn-.W';  . w-.-vi  r.  t hi  r-rnits  are  repca*,'-!  s,  : :w' 

a comjiarison  with  tlie  otlier  > -i.c  c --.i'  1-  . 

Viscous  lle::i  ll.  ' 

: . Ei'fect  of  Changi’.r  f i.*-  Fi  r ing  Fiti-f  io-:  Am;  , i t ude 

Tiie  magi.i  f ica’ ion  fa 'ter  varies  inversely  as  the  .-haf.ge  if. 
tile  forcing  functioti  jaim;- M ’ ,.j-- . 
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2.  Effect  of  Chariging  Mass 

Tlie  natural  frequency  and  the  damping  ratio  vai-y  inversely 
as  the  square  root  of  t.he  mass  change  . The  frequency  ratio  varies 
directly  as  the  mass  change  . 

3.  Effect  of  Chariging  Stiffness 

The  magnification  factor  and  natural  frequency  varies 
directly  as  the  stiffness  change.  The  damping  and  frequency  ratios 
vary  inversely  as  the  square  root  of  the  stiffness  change. 

Velocity-Squared  Damping 

i.  Effect  of  Changii  j the  Forcing  Furiction  Amplitude 

As  in  the  case  of  viscous  damping,  the  magnificat.ion  factoi- 
varies  Inversely  as  the  change  in  the  forcing  function  amplitude;  however, 
the  proposed  damping  ratio  varies  directly  as  the  change.  The  proof  of 
the  above  statements  is  demonsti-ated  by  solving  the  equation  of  motion  for 
F = 2.0  and  F = 1.0  and  comparing  the  results  at  several  different  fre- 
quencies. The  results  of  these  runs  are  show>i  in  Tables  7 and  8.  Table 
7 shows  the  results  of  the  run  with  F = 2.0  and  8 ^ - 0.001499089.  All 
parameters,  bot.h  steady- state  and  ti'ansietit,  ai-e  seen  to  be  exactly  equal 
to  the  resiuts  for  F = 1.0  and  8^  - 0.002998178  shown  in  Table  8.  It, 
would  be  absolutely  iricorrect  to  plot  the  response  diagram  with  the  vis- 
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the  change.  The  damping  ratio  varies  inversely  as  the  stiffness  change. 
The  proof  of  the  above  statements  is  demonstrated  by  solving  t.he  equation 
with  k = ■§■,  6^  = 0.002998179?  a forcing  frequency  of  7.071068  i-adians/ 

second  (Table  9)-  This  permits  the  result  to  be  compared  with  the  solu- 
tion for  Sj  = 0.005996358,  k = 1.0,  and  a frequency  oi'  10  radians/sec  na 
(Table  10 ).  The  following  changes  are  noted: 

a.  The  displacement  varies  inversely  as  Ak,  l.e..  1 

b.  The  velocity  varies  inversely  as  the  square  r < i f t!i'.' 
change , i . e . , /2 

c.  The  accelerations  are  equal 

d.  The  energy  var'es  inversely  as  the  change,  i.e..  . 

e.  The  phase  angles  at  the  positive  peaks  are  equal 

f.  The  time  varies  inversely  as  t.he  square  root  of  the 
change , i . e . , /’2  . 


These  changes  show  that,  the  coi’rect  expressi^ins  f r the  r'  S[  :.:•?• 
ratios  were  used. 

3.  Effect  of  Changing  Mass 

The  damping  ratio  varies  it.versely  as  the  :nass  change's  and  the 
frequency  ratio  varies  inversely  as  the  s<iuai'e  r c.t  1'  the  mass  ciiatige. 

The  effect  of  doubling  the  mass  f 't-  a forciiur  frequency  ‘f  7 1 Ot?-h 
radians/second  and  8^  - O.Oo,.  998179  is  sh  'Wn  irj  Table  II.  Tlie  rresj . .cu- 
ing data  for  a forcing  frequency  >9’  lo.U  raui  ans/sei'.  cid  and  9 . . v\  l-e<v\-.; 

is  shown  is  Table  Id.  A ■■•jm;  ari  s.'u  f these  tw  ■ figures  .'^h  w 

a.  The  di S{' 1 ac.-mcnt  ar-  <,-q  jal 

b.  The  vel  "'it.y  vari>  .'  inversely  a.-  tin  r.)  ,ai-  i'.  • t 1'  t.he 

clvange-,  i.e..  — ! — . 
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c.  The  acceleration  varies  inversely  as  the  change,  i.e.,  -g 

d.  The  energy  per  cycle  is  equal 

e.  The  positive  phase  angles  are  equal 

f.  The  time  varies  directly  as  the  square  root  of  the 

change , i.e.,  . 

D.  Adams-Moulton  Numerical  Analysis 

The  dimensionless  numerical  steady-state  resq'onse  cxrves  given 
in  Chapter  IV  were  plotted  from  the  data  gathered  by  the  Adams-Moulton 
method.  The  same  equations,  coefficients,  and  forcing  function  wei-e 
used  in  these  solutions  as  were  used  in  the  linearized  solutions.  The 
curves  are  plotted  using  the  dimensionless  ratios  given  in  Table  3- 

Since  the  concept  of  the  dimensional  analysis,  as  introduced,  is 
not  generally  used  in  the  analysis  of  dynamic  systems,  it  is  necessary 
to  give  a few  results  which  will  serve  for  elucidating  the  theoretical 
considerations  given  in  Chapter  IV.  The  results  are  limited  in  number 
because  of  the  computer  time  involved;  however,  these  data  are  suffi- 
ciently accurate  to  demonstrate  the  proposed  procedure. 

The  curves  are  generated  by  assigning  vai'ious  fixed  values  to 
the  damping  ratio  and  varying  the  frequency  from  5 to  ^0  radians/second . 
The  value  of  the  damping  ratios  covei-  the  range  above  critica'.  . i.e..  t.he 
point  where  the  aiiqWitude  decreases  cont  inuously  as  the  iiiput  frequency 
is  increased. 

As  in  the  viscous  case,  all  curves  show  a j rt.  .orainate  resonance 
peak;  howrver,  the  Coulomb  and  di spl acement- squared  -ases  inuicate  ivs,.- 
nances  in  additioi;  to  the  predomitjant.  peak.  Tlie  additiorial  } ■ aks  ai'” 
must  noticeable  in  the  fihase  angle  respoiise  curves  arn!  are  si'i-n  i 


n 

i, 

s 

dependent  on  the  darnping.  In  general,  for  small  values  , dampi.ng  a'd 
i 

low  frequency,  the  output  amplitude  is  substantially  the  same  as  t:>  -• 

At  the  higher  frequeircies  the  output  oscillations  decrease  asicr.j  t ’ ; -a. 
toward  zero. 

E . Time-Histor'ies  of  Motion 

Unlike  the  input  .motion,  the  output  motion  is  a dist  . r-teu  sirrss.  dra. 
function  of  time;  however,  bot.h  functions  have  the  sa;r,e  frequency  f.r-  a., 
damping  cases.  The  degree  of  distortion  varies  from  siriq  le  sinus..iaa;  - 
tortion  to  i-andom  distortion.  The  sinusoidal  dist,.i'ti.  n iras  a.mj;it  . 
frequency,  and  phase  components,  depending  on  the  type  and  de.yree  -.  f sa;:.: - 
ing.  In  general,  velocity- squai'ed  daiiipiiig  for  small  values  cf  the  ua;.'u  inr 
constant  is  frequency-distoi’ted. 

The  random  distortion  is  due  to  the  discontinuities  in  U,.-  uamj  In, 
function.  These  discontinuities  are  c,f  such  magnitude  that  it  is  v i :‘t  .a . . ,v 
impossible  to  attach  any  mathematical  descrlptiicri  t>;  the  irn  st  j na:.t 

disturbance. 

Although  one  is  intuitively  aware  oi’  the  ai 'rv'c  st  at>-;:ii  . i:  : ■ . 

not  Seem  Impossible  to  ti'ace  a dlsl.ort i. ,n  ^s  tlie  i'am[  ing  .a;;sta:d  !•  • 
increasingly  lar'ger  So  tiiat  a final  idei.’ if icat  iot;  ^f  the  chai’a  ; ■ 

would  be  iiossible.  A1  s.  . the  suci.’oss  '-p  ruci;  an  atia.ysis  woul.:  u'-;  ■ 
the  size  of  the  iricrement  al  stops  in  t.lu-  e>:pf  r i menl  al  piaco.i.ro.  A.'n  .,'i, 
sevei-al  -attempts  were  made  to'  t ak<  small  i'  ■romonl  s fr  ■!:;  a I:,  wn  • t.  : 1 • 1 ■ . 

the  variations  in  artipiitnio.  | ha.-’  . ati  i I’n-q'toncy  .hist  i • tt  ro  -j;,  : a;'*,  ■ 

di  scotit.  inui  ty  w>  *1  i r:  't,  a mean  ins  t‘ ai.-alysis.  I:,  ‘i.o  ■•r-i  ■a.'- 

t.he  amplit  in  . phaso.  an:  I'S'  j .■  • var;.  ’.or  a '.v  • o f tti  t i n:  ‘ •'  ''  r . 

the  t i:;io  , 1'  ■ irron.'o  f t.li..  Il:  ■ ; tiMijt;,-  ■ ;■.•!;  t i . a ’ -I 
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completely  different  set  of  initial  condiLions  occur.  Also,  it  is  i-m- 
possibie  to  ascertain  the  length  of  the  transients  induced  aftei-  the 

I 

discontinuity. 

The  time-history  curves  for  velocity- squarea  dampin,'.  Figs.  2.' 
through  2k  , have  the  same  general  shape.  The  shape  of  t.he  CoU-v-  was 
estimated  to  have  the  sa_>ne  general  characteristics  as  the  resp^;.se 
diagram,  i.e.,  for  leading  and  lagging  phase  between  the  inj'Ut  f.rcitig 
function  and  the  resulting  output,  the  variation  i!i  shaj..e  over  a -y-'ic 
of  motion  woulu  exhibit  a definite  frequency  characteristic.  A re; re- 
sentative  cycle  of  motion  was  chosen,  and  the  frequency  l\;i-  a quarter 
cycle  was  calculated  by  measuring  the  length  ol’  time  fr-,m  a ro  ■■r.;ssing 
to  a peak,  and  vice  versa.  The  phase  angle  vector  .liagram  is  shewt.  in 
Fig.  48  . The  {.>hase  angles  are  shown  i'.  r the  first  >iuarter  cycle.  1;.- 
tuitively  the  motion  should  vary  as  the  IVrcitig  functi'-ti:  therefer-  . a 
sine  wave  with  the  calculated  frequency  was  '-ssisnei.  Table  l4  gives  tia 
results  of  this  invest  igation  IVc  the  first  quarter  ef  la  tlcn  fr^;:. 


zero  to  the  peak  amplitude.  Tai  : e give: 
second  quarter  cy.tle  fr.^m  tlie  peak  > t:.e  . 
•uiations  for  the  ass  .;ned  :iu  tun  w-  :•••  : • rf 
The  values  of  the  n'mierh-a;  t i , n ai---  ■ 

close  cori't.-l  at  ii  the  vai  i" 
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Tf!E  AI/IPLirUDE  TIME  VALUES  OF  DISPLACEMENT  FOP  TliE  FIRST  QlViP.IER 
CYCLE  OF  MOTION  FOR  A SINE  WAVE  OF  30.0  RADlAI'iS/SECCND.  AN!  A 
PEAK  AMPLITUDE  OF  0.424875  INCHES 


X sir;  au;. 


uut 

Calc  ul  atfcrd 

0.075 

0.0318 

0.150 

0.063 

0.225 

0.094 

0.300 

0.  i:-x 

0.375 

0.  l.'M 

0.450 

0. 184~ 

0.525 

j.  .:M0 

0.600 

0. 

0.675 

0.  b5‘. 

0.750 

0. 

0.825 

0.  T - 1 

0.900 

0. 33uP 

0.975 

^.351? 

1.05 

U . 36  "'  > 

i.  125 

0. 3R34 

C 

0 

e.  3'^  'i' 

1.275 

O'  . ‘tCH 

1.35 

0.41L' 

1 .425 

0.  ’iU0'>4 

1 . 500 

0.4j3‘' 

1 . 575 

0. 4. 4v) 
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TABLE  14 


THE  AMPLITUDE  TIME  VALUES  OF  DISPLACEMENT  FOR  THE  FIRST  QUARTER 
CYCLE  OF  MOTION  ASSUMING  A SINE  WAVE  OF  28.690343872  RADIANS/ 
SECOND,  AND  A PEAK  AMPLITUDE  OF  0.42487';  INCHES 


Time , t 

(JUt 

X sin  (JUt 

Calculated  Values 

X sin  out 

Numerical  Solution 

0.0025 

0.07173 

0.030 

0.030 

0.0050 

0.14345 

0.0607 

0.0607 

0.0075 

0.21518 

0.090 

0.090 

0.0100 

0.28690 

0.120 

0. 120 

0.0125 

0.35863 

0. 149 

0.  149 

0.0150 

0.43036 

0.177 

0. 176 

0.0175 

0 . 50208 

0.204 

0.203 

0.0200 

0.57381 

0.231 

0.231 

0.0225 

0.64553 

0.255 

0.255 

0.0250 

0.71726 

0.279 

0..'79 

0.0275 

0.78898 

0.301 

0. 301 

0.0300 

0 . 86071 

0.322 

0.3-1 

0.0325 

0.93244 

0.341 

0 

Lo 

0 

0.0350 

1 .004l6 

0.358 

o.yi 

0.0375 

1 .07589 

0.373 

^ . i''- 

0.0400 

i. j476I 

0.387 

■ . ■ "* 

0.0425 

1.21934 

0. 398 

0.0450 

1 .29107 

0.4o8 

0. 41. 

0.0475 

1 . 36279 

0 . 4 1 5 

i .-1 

0.0500 

1.4345..' 

0.421 

' .4  1 

0.0';25 

1 . 506j4 

0.4. '4 

0.4.  ■ 
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TiiE  AMPI.ITUDE  TIME  VALUES  OF  DISPLA.CEMEIJT  FOR  THE  SEnOHD  QUARTER 
CYCLE  OF  MOTION  ASSUMING  A SINE  WAVE  OF  RAD  IANS/ SECOND 

AND  A PEAK  AMPLIT'JDE  OF  0.4ol+875  INC!!FS 


X sin  (ur  X siri  oi’. 


Time,  t 

Ujt 

Caicaiated  Values 

N'uuei'l  'ai  Si. 

0.0550 

1 . 72788 

0.420 

0. 4.0 

0.0575 

1 . 80642 

0.413 

V . 4 i H 

0.0600 

1 . 88496 

0.404 

0 . ‘f." 

0.0625 

1.96350 

0.39^' 

0.  3'  - 

0.0650 

2 . 04  _^04 

0.379 

0.  s8i 

0.0075 

2 . 12058 

0.363 

2.  4 

0 . 0700 

2.19912 

0.344 

. ■ U' 

0.0725 

2 . 27706 

0.323 

. ,'.■*  + 

0.0750 

2.35619 

0.300 

0.  -,01 

0.0775 

2.43470 

0.276 

t ■ . • 7' 

0.0800 

2.51327 

0.249 

o..;50 

0.0825 

CO 

C\J 

0 . 222 

0..  22 

0.0850 

2.67035 

0.193 

0.194 

0.0875 

2.74890 

0.163 

0. 

0.0900 

2.82743 

0,131 

0. 1 33 

0.0925 

c.%597 

0.099 

0.100 

0 . 0950 

. .98451 

0 

b 

0 . Oi  '* . 

0.0975 

3.06^305 

0.033 

0.03^ 

0. 1000 

3. i4l59 

0.0002 

0 . 000. 
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The  velocity  peak  is  sharpei’  than  the  viscous  case.  The  acceiei'atioii 
curve  has  a spike- shaped  peak  which  would  be  expected  for  the  skewed 
velocity. 

The  above  analysis  indicates  practically  no  incremental  jjhase 
variations  within  the  quarter  cycles  considered;  therefore,  the  differ- 
entiated displacement  should  generally  represent  the  velocity  and 
acceleration.  This  is  obviously  incorrect.  The  curves  and  the  data 
show  that  additional  terms  in  the  solution  are  requireu.  The  diffi- 
culties of  representing  arbitrarily  small  variations  in  the  displace- 
ment were  mentioned  in  the  statement  of  the  problem. 

These  variations  become  larger  as  the  forcing  function  magnitu.u- 
is  increased.  The  deviations  become  so  large  that  the  motion  becomes 
unstabLe  for  all  danqilng  cases,  since  the  response  is  a function  of  the 
forcing  mctgnitude.  The  unstable  character  of  the  response  may  occur  in 
the  displacement,  velocity,  acceleration,  or  in  the  higher  dei'ivatives . 

The  time-histories  for  Coulomb  and  displ  acement-squai-ed  damping  shew  th< 
effects  on  the  stability  of  the  motion.  Obviously  a similar  response  exists 
for  "jerk"  in  the  velocity-squared  case.  It  should  be  noted  that  the  dis- 
torted velocity  and  acceleration  have  tlie  same  period  as  the  di sp  1 aci-mei.t  . 

Motion  curves  that  are  typneal  of  each  respionse  curve  arc  giver, 
in  the  previous  chapter.  Both  the  steady-state  atid  transient  niot  ioiis  af 
correct  for  the  coi-responding  damping  case  and  I’an  be  used  for  the  de-- 
sign  or  analysis  ol'  any  system  desi-riled  1 y the  equati.ais  . f mot  i n.  The 
data  analysis  perfoj'med  by  trie  eii.mputer  | rogram  permits  the  j ritit.  ui 
tlie  pertinent  transient  informati'd  and  these  data  W'  n r n : f-  :•  caoh 

point  on  the  respons,e  diagram.’.  The  ini':  rmat  i.  n i .■  •,  ‘ pr'  Seii'ei  i.  -a 


of  tile  Volume  of  the  materi’il. 


CfiAFIER  VI 


COUCLUSIONS 

Linearizat  ioti  techniques  liave  been  applied  to  the  ci.ara  ■’ ••r:  sii  c 
damping  functions  to  allow  the  pj-ediction  of  the  steady-siate  response  cf 
the  systems.  The  amplitude,  velocity,  acceleration  and  phase  angle  diagrams 
are  compared  with  an  accurate  solution.  The  errors  for  these  assumptions 
are  given  in  Chapter  V.  The  magnitude  of  t.he  cj  rors  iniicate  that  ar.y 
value  of  the  damping  greater  than  zero  is  significant  and  .-ani-ot  bt 
"small." 

The  error  in  determining  a damping  coei'ficient  by  lineariz'  d 
assumptions  varies  over  a considerable  I'ange.  For  a fix'.i  fia  qu<-t; 'v,  * r.e 
error  approaches  zero;  however,  the  large  variations  in  the  w.av*'  fort,  ar  a 
function  of  input  forcing  function  magnitude  makes  it  impossible  to  give  a 
general  value  of  the  error. 

The  dimensionless  ratios  given  allow  one  to  r-ad  a'  a glance  from 
the  resonance  curves,  the  effects  of  the  resonant  fr' qta-nci.  r,  th-  ilrpla  a- 
ment,  velocity,  and  acceleration,  and  the  degree  of  the  lan,p;ng;  mcri„ver, 
the  motion  wave  forms  can  be  specified  directly.  Floltlng  * respctisc 
with  arbitrary  dimensional  quantities,  as  is  I'rcqu'cnt  ly  i.  ri'’  in  the  lit>  i - 
ature,  turns  out  to  be  less  survcyable. 

A difficulty  encountered  in  the  interpretatl on  of  d.vaiamic 
is  the  vagueness  of  the  effects  of  the  linearized  assumjd  ions . The  accu- 
racy of  these  linearized  assumptions  is  ('stabl  isiied  for  tdie  Kitz  ‘w  i-m 
and  equivalent  energy  methods  by  direct  cump'irison  with  the  numerii'al 
solution.  Variotts  limits  of  the  accuracy  hav  been  estaldirhcil  so  that 
the  design  problem  at.  hand  may  be  t rated  with  * lie  .lesir'\d  acctir-icy. 
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The  dimensional  analysis  of  the  problem  is  the  most  impori ant 
aspect  of  the  investigation.  Sufficient  I'a-miliariiy  with  the  dimetiSion- 
less  ratios  has  been  established  to  allow  a simple  and  convenient  pu'eset.- 
tation  of  the  results.  Simplicity  is  attained  by  looking  I'or  the  inde- 
pendent dimensionless  ratios  that  are  convenient  for  the  j resetit atlon  of 
universal  design  curves.  The  value  of  presenting  the  data  in  ter.ms  of 
the  viscous  critical  damping  ratio,  6,  is  definitely  incuj-rect.  The 
method  presented  can  be  established  before  determining  the  soluticti. 

This  method  of  establishing  valid  conclusions  from  a getieral  dimetisional 
analysis  provides  one  of  the  most  signiflcaiit  contributions. 

The  principle  advantage  of  studying  t.he  wave  f_:rmc  f ..  r the 
various  da.mping  cases  is  the  determination  of  the  predominant  freipuency 
content  in  tne  response.  This  will  per.mit  the  determination  of  iiai-motiics 
not  necessarily  in  the  environment.  With  t.he  frequency  range  est  alii  slice . 
exjieri-mentai  test  instrumentation  can  be  selected  with,  the  ji-.-per  response- 
characteristics . For  example,  the  j-roper  choice  of  a -ceJer.'mcter.  filtt-i-. 
amp]  i tier,  ati'l  recorder  response  couUl  not  be  made  t\:r  a iami  ■ d system 
without  some  knowledge  of  t.he  frequencies  involved.  Ai  s..; . a kn.  w 1 •.  ig- - .-f 
the  wave  form  is  invaluable  for  the  design  at.d  ai.a.ysls  of  --quipm- nt. 

The  steady- St. at and  transient  respo'nsc  -Irui'a -t-'i’i  st  ' -c  ..-f  th'-C-.- 
systems  with  other  for -irig  r.it.--tions  would  be  if  t rcni' ndo.-as  va;..e  ),  th- 

design  engineer.  A compl-  t.e  co  1 lect.i>.in  ^d’  su--!'.  .la' a exists  fo!-  vis  -. 

4 

damjied  systems:  how-'V-  r.  virtually  notiiing  is  t.ai  ..;i  ated  n---  r’dr.,'  l:;- 
response  of  the  nonlinear  sysiiins  givei.  iri  this  'hesis.  .'.oi-t  th'-  l‘  r-.-c 


function  is  completely  des<;ribed  iy  the 
the  dimension  1 >.'GS  ratios  c-m  le  a|pii-  1 


magti ; I u 1'  a;..i  sirq  . i'  t 

,i  1 ff.-t  l.v  t t i.-s.-  .-.vs*  - . 


1 

■ 

] 

I 


A 


I 

! ici 

future  work  might  include  investigations  of  combined  damping  cases.  The 
author  has  run  several  combined  damping  cases.  The  range  of  the  damping 
coefficient  can  be  predicted  by  the  Ritz  procedure  for  ail  cases.  In 
general,  the  response  does  not  contain  as  many  discontinuities  as  the 
basic  types  of  damping  presented.  It  should  be  noted  that  a general 
study  of  this  type  would  be  almost  impossible  because  of  the  number  of 
combinations  of  damping  coefficients  and  functions.  The  numerical  pro- 
gram in  the  appendix  can  be  used  for  a general  study. 
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